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Abstract. This paper quantifies the speed of convergence and higher- 
order asymptotics of fast diffusion dynamics on R" to the Barenblatt (self 
similar) solution. Degeneracies in the parabolicity of this equation are 
cured by re-expressing the dynamics on a manifold with a cylindrical end, 
called the cigar. The nonlinear evolution becomes differentiable in Holder 
spaces on the cigar. The linearization of the dynamics is given by the 
Laplace-Beltrami operator plus a transport term (which can be suppressed 
by introducing appropriate weights into the function space norm), plus a 
finite-depth potential well with a universal profile. In the limiting case of 
the (linear) heat equation, the depth diverges, the number of eigenstates 
increases without bound, and the continuous spectrum recedes to infinity. 
We provide a detailed study of the linear and nonlinear problems in Holder 
spaces on the cigar, including a sharp boundedness estimate for the semi- 
group, and use this as a tool to obtain sharp convergence results toward 
the Barenblatt solution, and higher order asymptotics. In finer convergence 
results (after modding out symmetries of the problem), a subtle interplay 
between convergence rates and tail behavior is revealed. The difficulties 
involved in choosing the right functional spaces in which to carry out the 
analysis can be interpreted as genuine features of the equation rather than 
mere annoying technicalities. 
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1. Introduction 



Long-time asymptotics of nonlinear diffusion processes have been a subject of 
much recent interest. The porous medium equation 

(1) p. = Iap- 

m 

is a prototypical example; it governs the evolution of a nonnegative density 
p{T,y) on [0,oo[ x R", as described in Vazquez book [IS] and its references. 
The basin of attraction IS] of its self-similar solution [IH] |1] [SH], and the 
rate of convergence of other solutions toittlljpJ][T9][36][2Q][3T][37][6][7] 
[9] [8] [22] have attracted a steady stream of attention since such rates were 
first obtained by Carrillo-Toscani [12], Otto [38] and del Pino-Dolbeault [T8] . 
Sharp results for convergence in entropy sense have since been extended to the 
full range of m G R by the quintet consisting of Blanchet, Bonforte, Dolbeault, 
Grillo and Vazquez [7], the quartet P] and the trio [^. However, few results 
are known concerning higher asymptotics, beyond the two improvements ac- 
cessible by choosing an appropriate translation in space (Carrillo, di Francesco, 
Kim, McCann, Slepcev and Toscani and the quartet [8] in various combina- 
tions and settings [13] [13] [31] [37] [10]) and in time by the Dolbeault-Toscani 
duo [22] (c.f. [3] [H]). In the one-dimensional porous medium regime m > 1, 
one has a full asymptotic expansion of Angenent [1] based on the spectral cal- 
culation of Zel'dovitch and Barenblatt [H]. Koch's habilitation thesis provides 
a potential framework for generalizing this to higher dimensions [32j. In the 
fast diffusion regime m < 1, a spectral calculation by Denzler and McCann 
[T^ [20] has been used to derive the first two corrections [2Z] [S] [22] to the 
leading order asymptotics [18] [3S] , but the higher-order modes never been suc- 
cessfully related to the nonlinear dynamics. It is framework for achieving such 
a relation which we develop for the first time below — inspired by ideas from 
dynamical systems, and sticking to the mass-preserving range m G ]^^, 1[- 

To achieve this, several obstacles must be confronted. The spectrum of Denzler 
and McCann [20] contains only finitely many eigenvalues below the continuum 
threshold, so only finitely many modes are in principle accessible. Moreover, 
there is an incompatibility between the spaces in which the eigenfunctions live, 
and the spaces in which the dynamics ([1]) turn out to depend different iably 
on their initial conditions. It is this differentiable dependency that we need to 
establish to justify the linearization which leads to the spectral problem. 

Guessing a space in which it will hold is far from trivial, however. One of the 
technical devices we use to achieve this appears independently in Bonforte, 
Grillo and Vazquez' [9] entropy-based approach to the special case m = 
namely, after linearizing the rescaled evolution in relative variables, we restore 
uniform parabolicity to this apparently degenerate equation by carrying out 
the analysis on R" viewed as a Riemannian manifold Ai with an asymptoti- 
cally cylindrical metric, known as the cigar. 
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Reconciling the above-mentioned incompatibility of spaces for the linearized 
and nonlinear dynamics complicates our analysis and yields a richness to the 
statements of our theorems which goes beyond the intrinsic complexity of a 
spectrum whose features include a web of eigenvalue crossings as m is varied. 
To obtain higher-order asymptotics, we need to work in weighted spaces which 
discount information appropriately at spatial infinity. The further we wish 
to penetrate into the spectrum, the more severe the required discounting. 
The corresponding linearized operators are no longer self-adjoint; they act on 
weighted Banach rather than Hilbert spaces. 

Since the qualitative behaviour of ([1]) varies considerable with m in different 
regimes, let us set 

(2) := 1 - 

n + p 

where the parameter p = — n is the moment-index introduced in [T^ . 
Thus m_2 = mo = ma = and m„ = for example. 

Herrero and Pierre [29j obtained remarkable local estimates for solutions to 
the fast diffusion equations if < m < 1 and mo < m < 1, which they and 
Dahlberg and Kenig [15j used to prove that every solution has an initial trace, 
and it can be extended for all t, and vice versa, there exists a global solution 
for every initial data that is a local measure. 

For m > mo, the evolution ([1]) of compactly supported initial data is mass- 
preserving; this is the regime we shall explore. For m < mo, the mass dwindles 
to zero in finite time; the basin of attraction and leading-order asymptotics 
describing this disappearance have been provided by Daskalopoulos-Sesum [16] 
and the quintet respectively. 

Although the dynamics ([1]) has no fixed point, a well-known rescaling using 
the self similar coordinates 

x = {l + 2pT)-^y, t = ^ ln(l + 2pT) , 

(3) /3 = (2-(l-m)n)-i = i(l + ^), 
u{t,x) = e("+p)"*p((e2p* - l)/(2p),e("+?')*a;) 

yields an alternate description 

(4) ^ = -Au"^ + -^V ■ (xu) 

CM m 1 — m 

of the same evolution, but with a fixed point in the new variables. This 
scaling coincides with the one used in our proceedings report pTOJ and by the 
quartet [S]. Conversely, we can express the density of the solution to the fast 
diffusion equation by 

(5) p(r, y) = {l + 2pT)-^'^u ln(l + 2pr), (1 + 2pT)-^y 
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The character of this rewriting depends on the sign of /3. If m > mo, then 
r — 7- oo is equivalent to t — ?■ oo and the r — )■ oo asymptotics for ([T]) translates to 
t — )■ oo asymptotics for (jl]). This remains true after a similar transformation 
for m = niQ. For m < mo solutions may extinguish in finite time and the 
asymptotics of ([1]) for t — )■ oo give information about the asymptotics of ([1]) 
at the extinguishing time, as in work of the quintet [7] and trio [9]. 

The stationary solution 

(6) Ub{x) = {B +\x\^y^ 
to (HD is related to the Barenblatt solution [15] |I] 

(7) Pnir,y) = i2pr + iy^^UB{{2pT + l)-''y), 

where B is determined by a quantity J podx called the mass, at least if this 
mass is bounded. It has been known to attract all solutions which share its 
mass since the work of Friedman and Kamin [25] ; see also Vazquez |13] . The 
Barenblatt solution has finite moments of exactly those orders smaller than p. 

It has been shown by Vazquez [ill Thm. 21.1] that weak assumptions like 



p{0,y)dy = j pB{0,y)dy 



and 

sup \p{0,y)/pB{0,y)\ < oo 



imply 

p{^,y) 



lim sup sup 

T— >oo y 
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pB[r, y) 

We investigate the sharp decay rate of the quantity | ■ | in (jH]), called the 
relative L°° distance. The known sharp rates of decay for integral expressions 
— such as \\p — Pb\\l^(r^) or the relative entropy \12\ [TB] [3S] [H] — imply a 
(non-sharp) rate of decay in relative L°° through the work of the quintet [7|. 

Because of the gradient structure with respect to the Wasserstein distance 
discovered by Otto [38] , convergence questions can be attacked using displace- 
ment convexity; see also [3Z]. It is within this framework that the linearized 
problem for m < 1 had been studied in great detail [20j. Unfortunately, there is 
no clean way of passing from the linearized operator to the full nonlinear equa- 
tion. However, by cleverly employing the entropy method, McCann-Slepcev 
[37] . the quartet [8], and the duo [22] were able improve on the sharp integral 
rates of convergence of Carrillo-Toscani (m > 1) [12j, Otto (m > m„) [38] 
and del-Pino-Dolbeault (m > m„) [IB], to extract first [37j [8J and second [22] 
corrections. The first order correction depends on centering the data, and the 
second-order on choosing a suitable translation in r. Our first result gives the 
sharp relative L°° rate of convergence for centered initial data. It improves on 
the rate found independently by the quintet [^ without centering, and implies 



the sharp entropy rate of convergence also found independently by the quartet 
for centered data [S]. 

Theorem 1.1 (Exact leading-order asymptotics in the relative norm). 
Fix < m G ]mo,l[ with mo = Suppose p{T,y) satisfies and the 

condition ([8]) holds for some B > 0. If m = 1712 = we assume in addition 
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(1 + \y\^)-'''^dy < 00 . 
Then there exists z G R" such that 
(10) limsup sup r 

Without the additional condition ([9]) in the case m = m2, we still get ffTOj) . 
except that the leading factor t is replaced by t^~^ , for any e > 0. 

Let us note that obtaining the sharp estimate (fTOj) without an e in all cases 
other than m = m2 is not automatic, but requires a detailed study of the 
linearized PDE, as given in the proof of Thm. 18.51 We rely on the manner in 
which the essential spectral radius depends on decay properties built in the 
space. The argument in case m = m2 and assuming the moment condition (Q 
is of a different nature, basically a spin-off of the theory and relying on 
self-adjointness in this case. 

The case m G ]mo, had already been proved by Kim-McCann [21], subject 
to a mild moment condition, which is not needed here, except for the similar 
moment-type condition iQ in the case m = m2. The case m > nin had been 
handled to order 0(l/r^^^) by McCann and Slepcev [37] in the somewhat 
weaker norm. Before this, Carrillo and Vazquez had given the sharp 1/r 
rate in the relative L°° norm for radially symmetric data [13], and a weaker 
0(l/r^/2) rate in the norm over the whole parameter range m > ruQ. The 
sharp rate for the case m2 < m < rrin remained open for a while, even though 
the role of the translation had been understood in the linearized setting by 
[20] and [13], before being resolved independently in the work of the quartet 
[8] and the present manuscript. The work of the quintet [7] gave sharp but 
implicit integral rates which extend even into the finite extinction range of 
parameters m < ntQ. The work of the trio [9j addresses the exceptional case 
m = m_2, in which the essential spectral gap f[T^ described below and found 
by Denzler and McCann [L9l ^20] vanishes. Whereas the works of the quintet 
and its subgroups are entropy-based, we use a dynamical systems approach 
that successfully bridges the functional analytic gap between the linearization 
argument of [20] and the nonlinear estimates, allowing higher modes to be 
accessed. 

The proof, which is found in Sec. involves preliminary results of fairly differ- 
ent flavors. A prominent role is played by a manifold A4 with one cylindrical 
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(9) 



Pb(0,?/) 



o„{t. v) 



< 00 



end, which is called the cigar in the two dimensional case; the same manifold 
was employed independently by the trio [5]. We express the equation in self 
similar variables and turn to the relative size u/ug as main dependent variable, 
as do the quintet [7j and their successors |H] [9] [22]. The equation for the quo- 
tient can be understood as a uniformly parabolic reaction/transport /nonlinear 
diffusion equation on the manifold M.. Well-posedness and smooth dependence 
on the initial data follow by suitable adaptations of known techniques. How- 
ever, since our setting deviates from more standard situations we provide com- 
plete proofs for the convenience of the reader, and for the convenience of having 
results tailored to our needs. The (rescaled and time translated) Barenblatt 
solutions are barriers, and local existence immediately implies global existence 
if m > mo = Below m < rriQ the time translation corresponds to an 

unstable mode, which is obvious by looking at the Barenblatt solution, where 
a time shift is used to adjust the r at which the solution extinguishes. 

The asymptotics of Theorem 11.11 are determined by the spectrum of the lin- 
earized operator. Its largest eigenvalue Aqo is zero. It corresponds to the 
rescaling 

Pa{r,x) = p{T/a'^,x/a) . 

The invariant manifold of this mode is given by the set of stationary solutions 
Ub, which is parametrized by B. For m > B is determined by the mass 
(L^ norm) of the initial data, which is a constant of motion. Once we adjust 
the mass, this spectral value becomes irrelevant. 

Equation ([T]) is spatially translation invariant. The corresponding eigenvalue 
of the linearization is Aio = —n — p. The invariant manifold is determined by 
all translations of Barenblatt solutions. For m > mi = the first moments 

^ n+l ' 

of the Barenblatt solutions are defined and conserved, and, by centering po we 
get rid of this mode. 

Equation is also invariant under time translations (which are equivalent to 
rescalings of space). The corresponding spectral value is Aqi = — 2p. This time 
there is no related conserved quantity, and there seems to be no direct way 
of determining the time translation parameter from the initial data |17] . This 
eigenvalue is responsible for the convergence rate given in Theorem 11.11 This 
suggests that we may improve the results of Theorem 11.11 by modding out the 
time shift as well, as was achieved in the entropy sense by the duo [22] , and as 
we hereafter independently achieve in the stronger senses of Remark 11.91 and 
Theorem 111.11 

Indeed, having quantified a rate of contraction, higher asymptotics become 
accessible, as we now explain. Accounting for the time shift is merely the 
first in a series of many improvements which are possible for m > m2, if 
we allow ourselves to measure convergence in weighted norms which suppress 
information in the tail region \y\ > and reveal higher modes. How strongly 

this information is suppressed depends on the degree of asymptotic accuracy 
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desired. It is necessary to introduce weights since Aqi turns out to coincide 
with the threshold of the essential spectrum in the unweighted Holder space 
C"(A^) of functions on the cigar defined by ( l33i) . For m < m2 on the other 
hand, any eigenvalues apart from Aqi turn out to be embedded in a spectral 
continuum and Theorem 110.11 shows we can obtain very rapid decay, but only 
for a more restricted class of initial data, since the weights in this case amplify 
the significance of tail information. 

Given p = 2(1 — m)'^ — n and non-negative integers k,i E'N, define 

(11) - Xek := {i + 2k)p + ni + 4:k{l ~ i - k) 
and 

(12) - Ar := (f + 1)' • 

In a certain critically weighted Holder space, the X^k for £ + 2k < 1 + p/2 
will turn out to be eigenvalues describing the exponential rate of contraction 
of u{t,x)/ug{x) towards the constant state 1, and Aq°'^* will turn out to be 
the threshold of the essential spectrum. For irrational m G ]m2,l[ and any 
A G ]Ao°°*,Aoi[, varying the weights in the linearization carried out below 
suggests there exist Ui^{x) depending on m(0, x) such that 



(13) 



{B + \x\^){u{t,x)/u,{x) - 1) - EA<.A<o%Ne-^* 



{B+ |a.|2)(p+2-V(P+2)'+4A)/4 



0(e 



C°'{M) 

as t — )• oo, where the sum is over integer-valued multi-indices l = {iik)i,k&i 
constrained so that zqo = = iio (assuming the mass and the center of mass of 
u{^,x) — Ub{x) vanish), and so that t-\:= J2Tk=o ^f-kX^k lies strictly between 
A and zero. The presence of essential spectrum also suggests that incorpo- 
ration of further terms u{x)e^^ into the sum cannot generally make the error 
term smaller than 0(e~^2"^^'* *). In this conjecture we see clearly the trade-off 
introduced by the choice of weight between weaker norms and the faster rates 
A of decay. This is a new feature to emerge from the present work, at least 
relative to the existing entropy-based results of the quartet |8] and duo [22j . 
and to the conjecture advanced in [20]. 

When m (or equivalently p) is rational, the coefficients Ui{x) would need to 
be replaced by a polynomial function uj^x, t) of time in case of an eigenvalue 
resonance, i.e. i; ■ A = ■ A with ^ i^k > 1, as in Angenent [1]. The possibility 
of such resonances can, for example, be ruled out either by taking m irrational 
(in view of the rational dependence of the eigenvalues f|TT]) on p), or by limiting 
ourselves to a A > 2Aoi within a factor of two of the spectral gap, which is 
Aqi = maXt i ■ A for m > 1712 and centered mass distributions. The restriction 
A > 2Aoi also simplifies the conjecture in two other ways, allowing us to prove 
the following theorem. First, it forces the functions uJ^x)/{B + \x\'^) to be 
eigenfunctions corresponding to the eigenvalue i; ■ A of the linearized operator 



fl23|) : this operator generates the long-time dynamics near the fixed point in 
the critically weighted Holder space C^^^{Ai) defined at flTUl) . Moreover, since 
the nonlinearity is analytic, its effects at this level are benign: the quadratic 
approximation of an orbit by its tangent is higher order than the accuracy 
A/Aoi < 2 accessible in the following theorem, one of our main results. 

Theorem 1.2 (Higher-order asymptotics in weighted Holder spaces). Fix p = 

2(1 - my^ - n > 2 (equivalently m e ]m2, l[) and A G [A™""*, Aqi] = [-(f + 
1)^, —2p] subject to the condition 2Aoi < A. Ifu(t,x) is a solution to (jlj) with 
center of mass and limt^oo \\u{t,x)/uB{x) — l\\L<x(^^n-^ = both vanishing, then 
there exist a sequence of polynomials {uek{x))ek, each element of which either 
vanishes or has degree £ + 2A; G ]1, | + 1[, such that 
(14) 



{B + \x\')iuit, x)/Ub{x) - 1) - EAKX^^KoM^y"' 



{B + \x\ 



(p+2-^(p+2)2+4A)/4 



= 0(e 



At\ 



as t ^ oo, where the sum is over non-negative integers /c, £ G N for which \ik 
defined by f llip lies in the interval ]A,0[, and for which i < 1 if n = 1. The 
functions Uik{x)/{B+\x\'^) lie in the \ik eigenspace of the linear operator ( l23l) 
on C"^^(A^), and the norm \\ ■ ||c«(A^) > II " ||l°°(R") is defined by fl33|) . 

Remark 1.3. The coefficients ueo{x) are given explicitly by Corollary \ 11. 31 
the invariant manifolds corresponding to these particular modes can be identi- 
fied using the same idea. 



Remark 1.4. The sum appearing in (]T4l) contains up to eight non-zero terms, 
depending on n, nip and on A. These correspond to eigenvalues from some 
subset of {Aoi, A02, A03, An, A12, A20, A21, A30}. 



Proof of Remark \1.4\ Formula f fTTj) implies the monotone dependence of Xik > 
Xik+i on k in the admissible range i -\- 2k G ]1, | + 1[, so X^k exceeds ]A,0[ 
unless —Xeo = i{p + n) < ip = — 2Aoi. This forces £ < 4. On the other 
hand £ + 2k < | + 1 implies {p + n)£ + 2A;(| + 1 - ^) < Xik- For this to be 

within 2Aio = -4p of the origin we need k < 4-£- < 4-£ 

which establishes the remark (except for the case n = 1 which needs to be 
argued separately). If n = 3 and m = rrin and A/Aio ~ 2 all eight terms may 
appear. □ 



These conclusions translate back to the original variables easily: 

Corollary 1.5. Fix p = 2{1 - m)-^ - n > 2 and 2Aoi < A G [A™"^*, Aqi]. // 
p{T,y) is a solution to with vanishing center of mass and satisfying ([8]), 
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then 



(15) 



>B(r,0) 
.pBiT,y] 



pB{T,y) 



E 

A<Affe<0 



tx,fc((l + 2pr)'^l/)) 
5(1 + 2pr)-^«/2p 



(p,(r, 0) /p,(r, y)) Vb+2)2+4A) /(2p+2n) ^ ^A/2p 



L°°(R") 

remains bounded as t ^ oo, the notation being the same as in Theorem \l.B . 
Proof. Evaluate (HM using 2pr = e^^* — 1 to get 



((i?+|a;|2)/5)(^+^-V(p+2)^+4A)/4 



0(r^/2p) 



as r — )■ oo. Setting a; = (1 + 2j9r) ^y, comparison with ([3])-(I7]) establishes 
(US]). □ 

Remark 1.6. From (|TT]) and /3 = |(1 + on/y products of powers of r^^^'^ 

and ofr^^^"^^ can arise as cof actors of the polynomials appearing in ( ITSi) . Since 
A/2p > —2, we can actually replace both occurrences ofl + 2pT by 2pT in ( |T5l) . 



Taking A = Aoi = —2p, the sum in ( |T5l) disappears and we recover an esti- 
mate in an unweighted norm, which is Theorem 11.11 By taking A < Aqi, we 
now concretely identify the leading order correction appearing in the weighted 
asymptotic expansion. This yields a corollary analogous in our framework to 
the entropic improvement found by Dolbeault and Toscani |22j . 




Figure 1. Hate 7 and weight 6 from Corollary \1.7\ and Theorem \10.1\ 
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Corollary 1.7 (Second order asymptotics modulo translations and dilations). 

Suppose 1 > m > 1712 = and p{t, y) satisfies the assumptions of Theo- 
rem and the center of mass of p(0, y) is at the origin. Then there exists 
Tq G R such that 

^,'PB{r,y)Y p{T,y) 



lim sup sup r"* , 

r^oo y \pB{r,0) 

where, for n > 2, we have 

_ [2-(n-2)(l-m)]^ 



1 



pB[r - To,y) 



< 00 



(16) 



and 



(17) 



7 
7 
7 



8p 

2(p-2) 

P 
n+p 
P 



8{l-m)[2-n(l-m)] 
4-2(n+2)(l-m) 



2-(l-m)n 
2 



2— (1— m)n 



if m2 < m < rriQ , 
if rriQ < m < rrin+A 
if mn+A < m < 1 , 



5 = ^Al-^) = W^-<^-^)) 



n+pV2 "I 4 

5 = = 1 - m 



^ = ^(|-l-v/(i-l)^-2n) 



z/ m2 < m < rriQ , 
^/ rriQ <m < rUn+i 
if mn+4 < m < 1 . 



n-l 



n+1 ^ 



For n = 1, the first case applies to m2 < m < rup^, and the third case to 
TTT'p, < m < 1, with = 2(v^ + 1), the middle case being omitted. 

Remark 1.8. On compact sets, or sets that grow at rate no faster that , we 
therefore get the improved convergence rate of 0{t~^), whereas an unweighted 
estimate that is uniform over all of R" cannot he obtained from this theorem. 



Remark 1.9 (Even higher asymptotics). Here 7 = A/Aqi and 5 correspond 
to the choice A = max{Ag°'^*, A02, A20} =: A in Corollary \1.5[ Corollary |i. 7| 
asserts that with this choice, a suitable translation of the solution p{t, y) in 
time makes the summation in (fT5|) vanish. Assuming p{t, y) itself denotes this 
translation and letting u{t,x) be the rescaled solution ([3]), hypothesis (11031) of 
Theorem \11.1\ is then satisfied, which allows us to access even more terms in 
the asymptotic expansion than provided by Theorem M.^ up to A > max{A + 

\ \ cont T 

Looking ahead at Theorem 18.1^ it will be noticed that m„ is the value where 
the eigenvalues Aio and Aqi cross, and that rUn+A is where the eigenvalues 
A20 = — 2]9 — 2n, All = — 3p — n + 4 and A02 = — 4p + 8 cross. 

Corollary 1.10 (Third order asymptotics and affinely self similar solutions). 

Fix m = rrip = 1 — 2/{p + n) with p > 4, n > 2, and let p{T,y) satisfy the 
assumptions of Theorem \l.l\ with the center of mass of p{0,y) at the origin. 
Then a traceless nxn symmetric matrix Sq and function cr(r) exist satisfying 
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{da / dry^"^ = cb det S(r), where T,{t) = Sq + o-{r)I > and cb is a constant 
depending on {n,p,B), such that p{t, y) = Us{T.-^/\T)y) detS-i/2(r) solves 
([1]) for large r and 



[18) limsupsupr 



pB[r,y) 
PB(r,0) 



< oo 



where ^ = -\u/2p and 5 = (f - 1 - + 1)^ + An) /{p + n). 

A more detailed interpretation of Theorem 11.21 will be given in the next section, 
after an exposition of the proofs with some more technical details. The full 
proof (relative to all preceding technical lemmas) will be given in Section [TTJ 
Let it suffice to say here that there is an intricate interplay between tail be- 
havior 6 and convergence rate 7 that is responsible for the occurrence of the 
weights in the theorem. 

A partner theorem for the case m G ]mo,m2[ is available below (Thm. 110. ip . 

Although the results contained herein fall short of a full-blown invariant mani- 
fold theory for the fast diffusion equation, they represent a significant advance 
in that direction. Invariant manifolds have a rich history of successful appli- 
cations in partial differential equations, as in e.g. Wayne |16] and the refer- 
ences there. However, much of that work has been devoted to semilinear heat 
equations, and cannot be directly adapted to the quasilinear equation now 
confronted. Wayne himself raises the question of whether it is possible to ex- 
tend such an analysis to nonlinearities which, though still semilinear, depend 
on derivatives. The present work develops the relevant ideas and techniques, 
and applies them to an example which may provide insights and a blueprint 
for the general problem. It complements the studies of the porous medium 
equation m > 1 by Angenent [1] [2] in one- dimension, and in higher dimen- 
sions by Koch [32] . Still there are some significant ways in which our analysis 
resembles that of Wayne ^6] or his results on the convergence of 2D Navier- 
Stokes dynamics to Oseen vortices with Gallay [28j: we study the size of the 
solution relative to the fixed point of an appropriately rescaled dynamics 
as they do (namely we work with {u — Ub)/ub) , and like them we use weighted 
spaces to shift the essential spectrum and reveal higher asymptotics in (fH]) . 
However, we are forced to work in Holder spaces (like Angenent [Ij, but on an 
unbounded domain), rather than the Hilbert spaces of Gallay- Wayne — or for 
that matter of McCann-Slepcev [37], the duo [22], or the quartet [^ — and 
with a spectral theory for non-self-adjoint operators, which is more involved 
than the linear problems analyzed in those works. Furthermore, there is a 
limit A G [Ag"'^*, 0] to the resolution which we can hope to attain from this lin- 
ear analysis, in sharp contradistinction to the complete asymptotics accessible 
in the problems addressed by Gallay and Wayne, and by Angenent. Except 
in the special case m = m_2 [S], we do not know how dynamical information 

beyond the threshold imposed by our continuous spectrum could be resolved. 
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2. Overview of Obstructions and Strategies, and Notation 



The basic idea to prove the asymptotic results is taken from dynamical sys- 
tems as outlined in [20], [16], or the references there: the eigenvalues of the 
linearization 'ought' to determine the rate of convergence to the equilibrium 
Ug. For example, when the analogous smooth finite-dimensional evolution 

= -V{x{t)) G R" 
is linearized around a fixed point V{xoo) = we get: 

- Xoo) = -DV{x^){x{t) - Xoo) + 0{x{t) - Xoo)^ ■ 

If a{DV{xoo)) = {Ai < ■ ■ ■ < A„} with eigenvectors (f)i, a naive expectation 
(albeit not entirely correct due to the possibility of resonances) is 

n n n 

x{t) = y: Q0.e-^'* + c^M^^'^'^' + E E E ■ ■ • • 

i=l i=l j=l i j k 

Notice however, that we may neglect the iterated summations and ignore res- 
onances if we are content with asymptotics to order 0(e~^^^*) as t — )■ oo. 
This is the strategy employed in our theorems. Differentiability of V{xq) or 
x(t) = X{t,xo) with respect to xq G R" is crucial. However for the infinite di- 
mensional dynamics of u, the choice of the right space becomes an issue. Even 
more than in the semilinear context [27], the issue is plagued by conflicting 
requirements from PDE theory and from functional analysis. 

The Hilbert space setting (a weighted Sobolev space W^p'^, which will be called 
W^'^ here) in which the spectral analysis was originally carried out in [20] is (at 
best) inconvenient to deal with the nonlinearity. Its geometric interpretation 
in terms of mass transport also restricts the parameter range for m to the case 
where Barenblatt has second moments. 

The singularity of the equation near is an obstruction to getting a smooth 
semigroup in many Banach spaces: While the parabolic evolution moves u 
pointwise away from the singularity immediately, it does not necessarily 
move it out of any norm neighborhood of 0, unless weights are built into 
the norm that are adjusted to the Barenblatt profile. On the other hand, as 
in Angenent [T], smoothness, rather than mere continuity, of the semifiow is 
crucial for the dynamical systems approach via spectrum of a linearization to 
have a significance. 

The 'relative' variable v = u/ug introduces the appropriate weight, and Vazquez' 
estimate ([H]) guarantees that after finite time, every solution is bounded away 
from in the norm of v. The singularity still shows up in the fact that the 
spatial differential operator fails to be uniformly elliptic in R", due to the de- 
caying Barenblatt whose power multiplies the Laplacian. But now R" can be 
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endowed with a conformally flat Riemannian metric, such that the Laplace- 
Beltrami operator with respect to this metric coincides with the linearized 
evolution operator in its principal part. Working consistently with respect to 
this metric restores uniform parabolicity of the equation. R" becomes a man- 
ifold Ai that can be illustrated by isometrically embedding it into R"^"*"^ like 
the surface of a cigar. 



cigar manifold 



The natural space in which to get a well-posed initial value problem for the 
nonlinear equation is a Holder space, and Holder quotients need to be referred 
to the geodesic distance on the cigar manifold. Likewise, for derivatives, unit 
vectors with respect to the Riemannian metric should be used to get appro- 
priately weighted higher Holder norms. It is in these C"(A^) spaces that we 
get smoothness of the evolution for the nonlinear equation. The detailed dis- 
cussion of the cigar manifold and the Holder spaces are found in Sections H] 
and[5l respectively. 

Section [6] develops some auxiliary heat estimates which are used in Section [7] 
to carry over the parabolic Schauder theory for linear and nonlinear equations 
from local results in fiat space to the uniform setting on the cigar. (Self- 
contained proofs for the estimates of Section [6] are deferred to an appendix.) 
The key point is that the cigar can be covered with coordinate charts in which 
the distortion between its Riemannian metric and a (local) fiat Euclidean 
metric is bounded with all derivatives, uniformly over all coordinate charts, 
even though the global distortion between the a priori Euclidean metric on R"- 
and the cigar metric is of course unbounded. Thereafter, the analytic local 
semifiow for the nonlinear equation follows just like the linear estimates by 
contraction mapping, exactly as it would be done in the textbook model cases 
on a bounded domain. Iterative usage of the local estimates can be controlled 
by a priori estimates via the maximum principle to get the global semifiow. 
Let it be stressed that the existence of a semifiow for the FDE has long been 
established by Herrero and Pierre [22] for Lj^^ data, our concern here is to get 
a flow map which depends differentiably on the initial data. 

Before comparing the nonlinear flow with its linearization, let us remember 
the obstruction given by Vazquez Thm. 1.3]: There cannot be a quantita- 
tive convergence rate to Barenblatt for general data. The idea behind this 
result is that one can start with small lumps of mass arbitrarily far apart, and 
the time until the spreading Barenblatts that arise from these lumps overlap 
significantly depends on their distance. So an infinity of lumps (with finite 
total mass) can produce arbitrarily slow convergence rates out of arbitrarily 
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large distances. The spectral analysis two of us developed in [20] is unaf- 
fected by this obstruction, because it is based on a linearization within Otto's 
formalism [38], which takes the distance between lumps into account in the 
form of transport costs. Of course, now we have to rephrase the spectral the- 
ory within the framework of the Holder space C°(A^), which will introduce 
some modifications. However, the key eigenvalues whose geometric meaning 
was understood from [20] are still in place, and so is their interpretation. We 
do have to cope with tail behavior issues, however, in particular for small m 
where the Barenblatt solution has few moments. The spectral analysis is done 
in Section [HI Much can be understood in terms of qualitative arguments of 
asymptotic analysis, only the precise values of the eigenvalues determining the 
convergence rates require the same explicit calculation as in [20] . We find it 
illustrative to get also the precise spectrum in various weighted Holder spaces 
out of the previously established formulas. 

Like for the norm, bounding the Holder space norm does not prevent us from 
starting with lumps that are far away from each other. We therefore cannot 
expect a global estimate of the form \\u{t) —Ub\\i < C(||'u(0)||2) exp[— at] even 
if we were to attempt to control the constant in terms of a rather strong norm 
||m(0)||2. In contrast to results of the quartet [8], our results are inherently 
local. It would be interesting to see if Wasserstein distance bounds could 
combine our spectral theoretic approach with more global information, but 
the functional analytic difficulties to make the 'formal infinite dimensional 
Riemannian manifold' approach work for the nonlinear problem have eluded 
us and appear to be pretty insurmountable. The proof of Thm. 11.11 is carried 
out as a rather immediate consequence of the flow properties and spectral 
analysis developed; see Sec. O 

Finally, we need to address the influence of tail behavior on the spectral theory. 
The Hilbert space in which our first spectral analysis was carried out carries 
a power of Ug as a weight, and allows a certain growth rate of functions in 
this space if m > m2, while enforcing decay if m < m2. The growth rates are 
powers of which is exponential in the geodesic distance s from the center 
on A^; and this affects the essential spectral radius. We introduce weighted 
Holder spaces C^, which allow a spatial growth \xY' (or equivalently cosh'' s 
in said geodesic distance from Sec[3|). The critical growth rate allowed in the 
Hilbert space is ?7cr = | — 1 where p is defined by ([2|). In terms of spectral 
theory, C°(A^) is therefore closest to the Hilbert space when m = m2, be- 
cause in this case rjcr = 0. The eigenfunctions of the linearization have their 
own characteristic growth rates (independent of m), which do not automati- 
cally coincide with the 'no growth' requirement that came from the need to 
keep the influence of the singularity m = of the nonlinear equation outside a 
ball around the Barenblatt in the Banach space. There is a trade-off between 
growth hypotheses and convergence rates, and once the obvious geometric in- 
variances (space and time translation) of the unrescaled fast diffusion equation 
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are modded out, finer asymptotics (like, e.g., from tlie effect that the diffusion 
brings anisotropic initial data closer to the isotropic distribution Ug) need to 
measured in norms whose weight adjusts to the tail behavior of the lineariza- 
tion in the corresponding direction in function space, a tail behavior different 
from the one of Barenblatt itself. This is what makes Thm. 11.21 and its corol- 
laries look somewhat technical. But there is good reason to believe that this 
displays a genuine phenomenon, not an artefact of the method. The good 
news is that, in fact, this discrepancy in tail behavior can be captured solely 
in terms of an adjusted weight in the norm in which the convergence rate is 
measured, but does not change the convergence rate itself. This fact arises as 
a consequence of the maximum principle for the linearized equation. 

Getting precise convergence rates requires an estimate for the semigroup of 
the form 0(exp[At]) rather than merely 0(exp[(A + B)t]). It is in this step 
that the tail behavior discrepancy (and its resolution by juggling the weights) 
become significant. This is the contents of Thm. 18.51 and its proof. The 
finer asymptotics provided by Thm. 11.21 and its corollaries are then proved 
in Section [TT] as a consequence in much the same way as Thm. 11.11 after a 
Lyapunov-Schmidt decomposition of the function space, with the exponential 
dichotomy property provided by the spectral gap. 

In the following, we will need to rewrite equation (jl]) in various forms, bringing 
to light the various aspects outlined above; different variable names help to 
distinguish the different versions and navigate between them. The notations 
collected in Table [1] will serve as a reference in this endeavour. 



3. The nonlinear and linear equations in cigar coordinates 

We approach the proof of Theorem 11.11 by an analysis in self-similar coordi- 
nates. 

To simplify the discussion we work with the relative density v := u/ug and 
calculate, using the PDEs for u and Ug, 

vt = u-^ut = ^u^^Aiu'^v"') + Y^Ug^V ■ (xugv) , 

M^iA(u>'") = M™~iA(?;™) + 2m-iVm^ ■ Vt;™ + m^1(Am™) v"" 

= u'^'^Aiv"') - 2^Vu'^'^ ■ Vv"" - ^u-^V ■ (xug) v"^ , 

B \ / 1— m B 1— m B V -H/ ' 

Ug^'V ■ (xUgV) = VUg^'V ■ (xUg) + X ■ Vf . 

Now we use Equation to get 



/ \ 2 1— nil i2 2 2 13 

Ug V ■ [XUg) = n — Ug \x\ = n — h 



1 — m 1 — m 1 — m B + \x 

15 



2 ' 



p solution of the unsealed FDE ([T]) - arguments t and y 

u solution of the rescaled FDE ^ - arguments t and x 

generic name for function on uniform manifold (Sec. [T]) 

Ub Barenblatt solution ([6]), with parameter B determining its mass 

ti =u/ub solves (pij) : generic function in heat equation estimates SeclB] 

V Variable of linearisation with respect to v about 1 

w w = V — I — u/ub — 1; generic Holder function in Sec. [5] 

Do, Wo initial data at time 

V tilde refer to conjugation with a power of cosh s 
u exception: u = uq + u in Sec. [7] 

vi, ui local functions on uniform manifold (pulled back to coordinate patch) vi = rjiui 

I index for coordinate patches (distinct from £!) 

rji partition of unity 

Xi coordinate maps 

rup value of m when Ub has moments below order p: rup = "^^^^ 

p moment parameter: n + p = 2/(1 — m) 

L linearization operator ((22|) 

conjugated linearization operator, see ([25l) 

L, L, Ijrj various operators with same principal part as L, 

Ai cigar manifold 

M generic uniform manifold (of which is a special case) 

r = \x\ 

s geodesic distance from origin: B^^'^ sinhs = r; Bcosh^ s = u^^^ 

C" Holder spaces (space or space-time) - Sec. |4] 

C" Holder space with weight allowing growth; see ([79| 

Cb space of bounded and continuous functions 

a Holder exponent 

/3 see dSl) 

Q spectral projection for L or L,, on a finite dimensional space 

P = 1-Q 

S (analytic) linear semigroup generated by L - analogous for S,, 

£ angular momentum quantum number - except for ([28| 

Affc , A™"* , A eigenvalue, onset of essential spectrum, spectral parameter 

i^ikjVik eigenfunctions in various settings: Thms. ISTT] and 15^ respectivelv 

uek = veku™~^ 

Yip spherical harmonics 

ri^ricr parameter for conjugation of operator (growth in space); rjcr = 2~ ^ 

b°°{ri) coefficient of 1st order term of L,, at cx) 

c°° (77) coefficient of 0th order term of at 00 

T finite time step; we consider time-T maps 

T related to spectral parameter A, see Fig. [2] on page [51] 

Table 1 . Overview of notations used 

and consequently 
(19) 
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This reaction / transport / nonlinear diffusion equation will play a central role 
in our analysis. It will turn out that it is a uniformly parabolic problem (for 
< V < c) on an unbounded manifold with a cylindrical end, known as the 
cigar. Its structure can be seen more easily in polar coordinates. We may and 
do normalize the parameter B to B = 1. 

With \x\ =: r, we can write A = d"^ + ^^dr + ^As, where As is the Laplace- 
Beltrami operator on the unit sphere. So we want to rewrite (fT9l) as 



(20) 



n-l _L ^ + 1 )AsW"" 



m \ J \ ar^ r or I m \r 

m dr ^ dr^ ] ~^ m^r^ "^"^^ As^ 



[V — V 



2(l+"i)^ I n-1 



l—m 



r + 



in — l)r 



or 



J 2_^_9_^ 

l—m dr 



^ \l-mJ 



n(l— m) 



l + (l + r2)-ij(w-w"^) 



We set r = sinhs, hence 4- = jf^ = (coshs) ^ ;|- and 

' or or as ^ ' as 

[(l + r2)|:^i = ^(coshs)-if ((l + sinh^)(cosh.)-i ^t;™) 
= -^v"" + ^ tanh s ^v"" . 

m OS'' m as 

This s will be the geodesic distance from the origin in the Riemannian metric 
of the cigar manifold introduced below. We can rewrite equation ( !20l) as 



(21) 



Vt 



:V"' + - 



2(n-l) d 



m \ds'^ sinh(2s) ds 



v^'+itanhsy^Asv'' 



+ —[n 
m 



2(m + l)\ , d ^ 

— tanhs— 

1 — m / as 

2 , , ^ 5 4 

+ ^ (tanh s)—v + 

l — m OS (1 — my 



l — m 



n- 



1 + (cosh s) 



,v — V 



It is convenient to introduce the operator L as the linearization about 1 of the 
operator in fl2T|) . Namely, 



(22) 



L;=52 + 2(n-l), , 



sinh 2 s 
2m 

+ 1^- ^ 



ds + (tanh s)"^ As 



tanh s ds + 2n tanh s . 

m/ l — m 
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For later reference, we also calculate L in cartesian coordinates from f|T9|) : 

m ia 2(1 - 2m) „ / 1 „ 2 , ,2\„ 
(23) s 1-m V ''1-m''/ 
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Using L, we will also find it convenient to rewrite (1191) for w = f — 1 as 

2 2 / 2^^^™ \ 

(24) = L/i(w)H a;-V(w-/i(t(7))H U — \x\^]{w~h{w)) , 

1 — m 1 — mV 1 — m / 

where /;.(ti') = '■''''''^ = w + O(w^). 

Later on we will conjugate operator fl2^ by cosh'' s. In order to keep these 
calculations at one place we calculate 

cosh"'' s o ds o cosh'' s = ds + r] tanh s 

and 

cosh"'' s o dg o cosh'' s = 9^ + 2// tanh sdg + r]{r] — 1) tanh^ s + r] , 
which gives after some calculation 

:= cosh"'' s o L o cosh'' s 

5.^ + ^^^9. + (tanh.)-2As 
smh Is 

(25) _2(^-f-l-r7)tanhs9, 
+ (T^-f-l-^)'-(T^-f + 1) 

Since the leading order terms coincide with the Laplace-Beltrami operator 
fl29|) . we can see — as a Schrodinger operator for a quantum particle moving 
in a potential-well on the cigar manifold, perturbed by a transport term which 
spoils self-adjointness. The potential-well has a universal (coshs)"^ profile of 
depth 

d{r]) = (2(1 - m)-^ - r/)(2 + r/) = (p + n - r/)(2 + r?) , 
and is asymptotic to the constant 

(26) - c~(r/) = (2(1 - m)-^ - n - ri){2 + 7]) = {p - ri){2 + r]) 

at the s — )■ cxD end of the cigar. The transport term shifts mass along the cigar 
with outward velocity asymptotic to 

(27) - b'^{r]) = 2((1 - m)~^ - r] - 1 - n/2) = p - 2r] - 2 

as s — )■ oo, making it harder for the operator to support eigenstates unless 
r/>|-l. 
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The value rjcr = ~ 1 ~ 2 ^ ^ plays a special role, causing the transport 
term to disappear and restoring the symmetry property 

/•oo /» roo I* 

I I (L^^^-Ui) tanh'^~S(iw„_ic/s = / / -Ui (L^„.'i;2) tanh"~^ s (iw„_i(is 
Jo J§"-i Jo J§"-i 

of the operator L.^^^. with respect to the volume element fl5U]) on the cigar man- 
ifold. Equivalently, Lq is symmetric with respect to m™ times the Euclidean 
volume element in R", which is the analog of the self-adjointness in [20]. 

This same value r/cr = ~" t ~ 2 ■'■ Kiininiizes c°°{7]), producing a quantity 

which gives the onset of the continuous spectrum, the negative of what was 
called Aq""^* in [20], see Thm. 18.11 below. It combines with depth d{r]cr) to 
produce a vanishing ground state energy for the Schrodinger operator L^^^ 
on L^iM); c.f. Thm. 1121^1). We note that r/cr = iff m = ms = it is 
positive for m > m2 and negative for m < m2. This is the reason for the case 
distinction between Thm. 11.21 (proved in Sec. [TTj) and Thm. 110.11 

In the remaining sections of this paper we will explore the consequences of our 
calculations. We shall see next section that the metric defined by the leading 
part of our operators has a very natural interpretation as a Riemannian metric 
on a well-known manifold, the cigar manifold. 

The nonlinear differential equation f l2T|) is (for v bounded above, and bounded 
away from 0) a uniformly parabolic equation on the cigar manifold, and well- 
posedness is a standard consequence of estimates for the linearized equations, 
which is basically standard parabolic regularity theory. Therefore we expect 
that the asymptotic behaviour is controlled by the spectral properties of the 
linearized differential equation. But since the cigar manifold is not compact, 
there is nontrivial essential spectrum for the linear semigroups S{t) = exptL 
and Sr,(t) = exptL^ defined through ( 122|) and ( p5|) . One may easily read off 
the essential spectral radius: It is e^°°^^^^ and it depends on 77. 

We are mainly interested in modes related to eigenvalues outside the essential 
spectrum. All eigenvalues and eigenfunctions have been determined by [20] . 
and we can use these results, with due adaptation, which is done in Sec. El 
In case r] = rjcr, we can already anticipate that the number of eigenvalues 
will be finite and increase with the depth d{r]cr) = (| + 1 + ^)(| + 1) of 
the potential-well. In one-dimension the Schrodinger operator L^^^ has been 
studied in connection with transparent scattering, and k{k + l)(coshs)~^ for 
/c e N are known as the Bargmann potentials in that context; a reference for 
these is Sec. 2.6, Exercise 11 of [35], as well as Sees. 2.5 and 3.5 there. In our 
case = |. As mentioned in Section [21 we will need the conjugated operator 
to accommodate the growth of relevant eigenfunctions. The dependence of the 
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essential spectral radius on rj will turn out to be useful in estimating sharp 
decay for the linear semigroup. 



4. The cigar as a Riemannian manifold 

The cigar is an analytic Riemannian manifold which we denote by A^. It can 
be described as R" equipped with the metric 

(28) de%i := (^l + \x\^^ ^^dx'^i = {ds^ + tanh^ s dll) , 

where di'^ is the length element on the unit sphere and sinhs = \x\ was 
already introduced above. It is immediate from this formula that s is the 
geodesic distance from the origin. 

In two dimensions, the cigar is a soliton solution for the Ricci flow (see ch.2 
of [Ej). For the fast diffusion problem, the cigar metric is appropriate for 
various reasons: The leading part of L (l22l) coincides with the Laplace-Beltrami 
operator 



(29) 



. n/2 / X -n/2+1 

od,o(l + \x\^] odi 



5.^ + ^^^9. + (tanh.)-2As 
smh 2 s 



which is also the leading part of the operator Yl^=i-^i^ where Xj are the 
obvious orthonormal vector fields 

/ \ 1/2 

Xi:= [l + \x\^] Bi 



and Si := d/dxi is the standard basis. Specifically, 

n 

^ Xf = Am + {n-l) tanh s 



i=l 



As a consequence, we shall see that the parabolic estimates become uniform 
when distances are measured with respect to the cigar metric. 

The Riemannian volume element is 



(30) dfi := tanh" ^ s dsdun-i ■ 

The vector fields Xj do not commute, but their commutators of any order can 
be written as a bounded linear combination of these vector fields. For instance, 

[^-^^■] = (l + |^|2)l/2^i - (l + |^|2)l/2^* • 
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5. Uniform manifolds and Holder spaces 

We will need uniform Scliauder estimates for parabolic equations on Ai. The 
possibility of such uniform estimates relies on the fact that locally, Ai can be 
mapped into R" with bounded distortion, the bound being global, even though 
a global map with bounded distortion is not possible between and R". It 
seems expedient to elaborate on this principle in its natural generality. So we 
will prove these estimates for parabolic equations on a uniform manifold M, of 
which the cigar M. is the example we are interested in: We study manifolds M 
(not necessarily Riemannian) with a distance d which turns M into a geodesic 
space: A metric space M is called a geodesic space if, given two points a;, y G M 
there exists a path 7 : [0, 1] — )■ M from x to y such that (i(7(s), 7(t)) = |s — 
t\d{x,y). This requirement is in particular satisfied for the geodesic distance 
of a connected closed Riemannian manifold. 

Definition 5.1 (Uniform manifold). Let M 6e a manifold with a metric d 
which turns M into a geodesic space. We say (M, d) has a uniform struc- 
ture (resp. uniform analytic structure) if there exist two constants R > and 
C > and coordinate maps Xx '■ Bfi{x) — ?■ R" for each x G M such that 

(31) C-H{y,z)<\xM-Xx{z)\<Cd{y,z) for all x e M, y, z e Br{x) 
with (resp analytic) coordinate changes Xy ° Xx^ satisfying 

(32) \d^{xy o Xx')\ < ^n Xx{Br{x) n BM) 

for any multi-index /3 of length < k (resp. for any multi-index) . We call R 
(and C) a radius (and constant) of uniformity for M. 

The first condition says that balls of radius Br are uniformly bilipschitz equiv- 
alent to subsets of R*^. The second condition (15^ implies that the coordinate 
maps are uniformly (resp. analytic). 

For our purposes, would sufficient regularity, but the cigar manifold is even 
analytic; and it is not hard to see that the exponential map provides such 
coordinate maps for the cigar manifold with R = 1 and a suitable constant C. 
Observe that we do not require a Riemannian structure on the manifold. But 
if the manifold carries a Riemannian metric then we use the associated metric, 
which turns the manifold into a geodesic space. 

The following lemma implies that the volume of balls is at most exponential in 
the radius. Given an open subset A of M and r > 0, we denote by N{A, r) the 
maximal number of disjoint balls of radius r in A (if such a maximum exists; 
else N{A,r) := 00). 

Lemma 5.2 (Packing r-balls into p-balls). Suppose that ISA is a uniform man- 
ifold of dimension n with a radius and constant of uniformity R and C . Then 
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the following bound holds 

iV(BA.).r) < I e?")" , >Jr<P<R, 

[ {max{C^R/r, 1})"(5C2)5"p/-R if max{r, R} < p .. 

Proof. Suppose first that r < p < R and consider a ball Br{xi) contained in 
Bp{xq) C Br{xq). Then 

Comparing the volumes of the images we see that there can be at most 
(C^p/r)" such disjoint balls in Bp{xo). 

Secondly, let's suppose r — ^, and p arbitrary. The previous case p < R 
will serve as the beginning of an induction with steps from p to p + r. So 
suppose {Briyj) | j = 1, . . . ,iV} is a maximal packing of r-balls in Bp^xo). 
Then the B2r{yj) will cover Bp_r{xQ) as a consequence of the maximality of 
the packing. But then the B^riVj) will cover Bp+rixo)- This is because for any 
z G Bp^r{xo), we can find a z' E Bp^r{xo) that has distance 2r from z; and 
then z G B/iriUj) whenever z' G B2r{yj)- Having constructed this covering, 
we now take a family of disjoint balls {Br{zi) | i G /} in Bp^rixo). As each 
Zi must be in some B^riVj), the Br{zi) lies in B^riVj)- As there are at most 
(5C^)" such Br{zi) in each given B^riUj), and there are only N many y^'s, we 
conclude N{Bp+r{xo),r) < {5C^)''N{Bp{xo),r). 

We have thus proved inductively that N{Bp{xo), R/5) < (SC^^^np/R 

Now we turn to the general case: If r > the estimate for r = is still valid 
trivially. If r < suppose we have a maximal collection of Br{xi) in Bpixo). 
We also place a maximal collection of Bji/r,(^yj) in Bp{xo). As before, the 
maximality then ensures that the collection of B^ji/^dij) covers all of Bp(xo) 
and each Bj.{xi) is contained in some Br+sR/^d/j) C B4^jij^{yj). 

By our second estimate, there are no more than (5C^)^"^''/^ of the BRi^{yj) 
fitting in Bp[xo), thus bounding the number of yj, and by the first estimate 
no more than [C'^R/rY of the Br{xi) fitting in each BiRi^{yj). □ 

As a consequence, given r < R/3 there is a sequence of points {xj) with 
distance at least r such that M is covered by the balls B^rixj). There are 
coordinate maps Xx^ defined on these balls since r < R/3. No point lies in 
more than (3C^)" of the balls B^j-ixj)- We fix r and the points Xj in the sequel. 
There is a partition of unity {rjf) subordinate to this covering with uniform 
bounds on derivatives (up to fixed order) of rji o Xx^ on Xxi{Br[xi))] we refer 
to this property as 'uniform smoothness'. 

We may choose such a uniformly smooth partition of unity. The structure of 
uniform manifolds allows to construct useful functions, specifically a smooth 
approximation to the radial coordinate (i(-,a;o): 
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Lemma 5.3 (Approximate radial coordinate on a uniform manifold). Choose 
Xo G M. There exists a function p : M — )■ R with bounded derivatives such 
that p{x) — d{x, xq) is bounded. 

Proof. We define the functions pi on the l^^ coordinate chart Xxi{Bj^{xi)) as 
regularizations of {rji o x^^) ■ {d{Xx^{-),Xo) ~ d{xi,Xo)), and then we can take 
P = Y.Vi-ipi°Xxi+d{xi,xo)). □ 

Given a uniformly smooth partition of unity rjf (as explained before the lemma) 
on the uniform manifold, we define Holder norms 



C«(M) 



snp\\{riif) o X-x^Wc-ix.ABRixO)) 



where on subsets X C R", 

(33) ||/||c;.(x):=max{[/],,||/|U^} 
with 

0<|x'-X'|<l F 

This definition depends on the choice of the partition of unity, but different 
choices lead to equivalent norms. The presence of the L°°-norm in (133|1 permits 
to drop the constraint |x' — x| < 1 altogether in the definition of [/]„, or else 
to replace it with a different bound |x' — x| < R, again providing equivalent 
norms. Likewise, it is easy to see that an equivalent norm can be defined 
without reference to a partition of unity by 

(34) ll/llc<^(M):=max{||/|U^(M),[/]°} 
with 

[/]"•= Slip ^ , ■ 

0<d{x',x) "'{X , X) 

We will use the variant f l33p for the general proof of the Schauder estimates, 
but the equivalent variant ( l34l) . when dealing with Ai specifically. In this case, 
d obviously refers to the cigar metric on Ai, not the euclidean metric on R". 

We can similarly define the spaces C^'"^ by the maximum of the norms (133|) 
of all the derivatives up to order k. 

Holder spaces C°([0, oo) x M) on space-time cylinders can be understood using 
the parabolic metric 

dp{(ti,xi) ; {t2,X2)) = max{|t2 - ti\^ ,d{xi,X2)} ■ 

This leads to 

||/||c-([o,r]xM) := sup \\{r]if) o X~lc"([o,T]xx.;(Bfl(x,))) > 
where on subsets X C R" 



(35) 



||/||c-([o,T]xX) := max{[/]^;„, [f]t;a/2, ||/||l-} 
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with 



[f]x;a '■— SUPt SUPo<|i:'_a;|< 
[/]t;a/2 := SUp^. SUPo<|t/_t| 




Similar comments about equivalent norms in the style of apply in an 
obvious manner. In the literature, these spaces are sometimes denoted as 
C"/^'"", and slight differences in the definitions lead to different, but trivially 
equivalent, norms. 

We will henceforth find the abbreviations 



(36) Mt := [0, T] X M , Mr := [0, T] x M , := [0, T] x 



Using multi-index notation /3 := (/3o, • • • , /Sn) with the parabolic weight 



\(3\ := 2/3o + /3i + . . . + /3„ and := Sf^xf^ ■ ■ -Xf", we let, for functions 



and can obtain C^'" norms on My or on [0, oo[ x M by partitions of unity as 
before. 

We note that not all authors include the term [d^ f]t-(i+a)/2 in their definition 
of the parabolic Holder norm (e.g., Krylov [33] and Friedman [25] don't, but 
Ladyzhenskaya et al. [31] do). For k = 2, this term is indeed controlled by the 
other terms, as can be seen in Krylov's book from his Ex. 8.8.6 in connection 
with the proof of Thm. 8.8.1. The argument will generalize for even k. The 
non-equivalence of the two styles of parabolic Holder norms for odd k will not 
be an issue for our purposes. We find it convenient to repeat the simple norm 
equivalence argument here, for easy reference. It suffices to do the argument, 
which is local, in flat space because of the local equivalence of the metrics. 
For the rest of the section, we simplify notation by giving the arguments for 
one space dimension (writing a scalar x), with the generalization to higher 
dimensions being obvious. We estimate 

\Wr,(t,x) - Wcc{s,x)\ < 



useful. 



on R", 



ll/llc*,- := max <^ sup ||9'^/||c-, sup sup [<9^/]t;(i+a)/2 




< Wa:(t,x) 




w{s, X + h) — w{s, x) 



h 



+ 



w{t,x + h) — w(t,x) w{s,x + h) — w{s,x) 



h h 




+ \t - s 



wt{9,x + h) - wt{9,x) 
h 
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where in the second term the mean value theorem was apphed to the function 
g{t) := "'(*'^+^^~'^(*'^) _ Using the Holder estimates for Wxx and Wt, we conclude 

\Wx{t,x) - Wx{s,x)\ < h\t - [w^^]t-a/2 + \t - [Wt]x;a ■ 

We now obtain \wx{t,x)-Wx{s,x)\ < \t-s\^^+°'^^'^{[wxx]t;a/2 + [wt]x;a) by letting 
h:= |t-s|V2. 

We will heavily rely on spaces of Holder continuous functions. They have nice 
algebraic properties: products and composition have good properties in these 
spaces and there are optimal estimates for linear parabolic equations in these 
function spaces. We begin with the algebraic side, for use in the next section. 

Lemma 5.4 (Estimates for some nonlinearities in Holder spaces). For D C 

open and convex, h G C^"*"^(D) and a pair of functions {f,g) with range in a 
compact subset of D, the following estimates hold on M. and on M.^: 

(37) \\ho{f,g)\\c..^<C{\\f\\cK.A\g\\c^..) 

where hounds on h and its derivatives on the convex hull of the range of f and 
g enter into the dependence of C on the norms \\f\\ and \\g\\; 

(38) ll/^llc^." <c||/|U^||^?||c^..+c||^7|U^||/||c... ; 

(39) \\h o (/,^)/||cfc.- < C(||/||cfe,a, ||^||c*.«)||^||L-lkllc*." • 
// h is analytic, then the map 

(40) C''' X C'-'^ 3 if, g) ^ h{f, g)g' G C"='" 
is an analytic map of Banach spaces. 

Remark 5.5. The lemma will he applied in a slightly more general situation 
with h depending also on x and t. This however follows from the lemma ahove 
hecause the statement remains correct if we consider vector valued functions 
f and g, and we may replace f hy the vector {t, x, f) for hounded functions x 
and t. Typically we will consider local coordinates with uniform estimates for 
all quantities which occur. 

Remark 5.6 (Notions of Analyticity). There are different equivalent notions 
of analyticity of maps defined on open sets in real or complex Banach spaces. 
We refer to section 15.1 in chapter 4 of [H] for a detailed exposition of the 
facts stated helow. Let X and Y he complex Banach spaces and U C X an open 
subset. A map f : U Y is called weakly holomorphic if it is locally hounded 
and if for every one dimensional affine subspace L C X and every I G Y* the 
map L nU 3 X i — y l{f{x)) G C is complex differentiahle. It is an immediate 
consequence (see exercise 2 in section 15 of\n\) of the Cauchy integral formula 
that a weakly holomorphic function is continuously differentiahle and complex 
differentiahle. Locally hounded complex differentiahle functions are analytic in 
the sense that the Taylor expansion converges in a hall. Real analyticity can 
he defined in terms of convergent Taylor expansions. Real analytic functions 
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can be extended to complex analytic functions, and conversely, any restriction 
to a real subspace (if X is the complexification of a real Banach space) of the 
real part of a holomorphic function ( assuming that Y is the complexification of 
a real Banach space) is analytic. Moreover both the real analytic and complex 
inverse function theorem holds. In the sequel we do not distinguish between 
different definitions of analyticity and use whatever is convenient. The results 
of this paper do not depend heavily on the notion of analyticity. Its main pur- 
pose is a nontechnical approach to differentiability and higher differentiability 
in the previous lemma. 

Proof. The first inequality fl37j) requires a standard calculation. We may re- 
strict ourselves to h o f and verify 

\\hof\\c.,.<C{\\f\\ 

by allowing vector valued functions /. It is a routine induction over that 
for any multiindex /3 with < /c all terms of d^h{f) can be estimated in the 
desired way. 

We get the k = case of ( l38l) (which is the most important for us) from 

f{x)g{x) - f{x')g{x') = f{x){g{x) - g{x')) + {f{x) - f{x'))g{x') . 

The case A; > 1 of (155]) is more subtle and is of an interpolation type, since 
naturally occurring intermediate norms as, e.g., ||5f||(70,Q must implicitly 

be controlled by the extreme norms on the right hand side. It suffices to prove 
f p8|) for fixed compact support, since the general case can be reduced to this 
special one by means of a locally finite partition of unity. In that case however, 
due to the local equivalence of metrics, one can prove the estimate in fiat 
space. A tedious, but elementary, proof can be given inductively. For instance, 
to estimate a term like [fi'la, one can use interpolation inequalities as 

Krylov's [33] Thm. 3.2.1 in the elliptic case, or Thm. 8.8.1 in the parabolic case; 
the scaling weight e contained therein will be chosen in terms of the norms, 
such as to remove one higher norm from products like E ^"||/||c'i'"||5'||c^'"- 

A less tedious, but not as elementary, proof in the case of R" can be found in 
the book of Runst and Sickel |4]J, Sec. 4.6.4, and the argument there gener- 
alizes to our setting. (Note that their setting in Besov spaces includes Holder 
spaces C'^'" = B^^.) The Littlewood-Paley decomposition (dyadic partition 
of unity in Fourier space) that enters implicitly already in the definition of 
Besov spaces (pg. 8 of [41j) needs to be made anisotropic to refiect the par- 
abolic scaling. We omit tedious details, which would only distract from the 
ideas of the present paper. 

Inequality ( 139|) is a consequence of the previous bounds: 

11^ o (/,fl')fi'lcfc>- < c||/io (/,fi')||L-||fi'||L-||5'llc"=.- +c\\ho (/,5')||c'=.-llfl'llL- • 
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Now suppose that h is analytic. It can be extended to a holomorphic complex 
map. We complexify the Holder spaces by considering functions with complex 
values. Left composition with h is clearly a weakly holomorphic map, hence 
holomorphic. The restriction to real valued functions is therefore real analytic. 

□ 



Remark 5.7. The reason why we insist on fl38|) rather than the simpler al- 
gebra property \\fg\\ < c||/|| \\g\\ in C*^'" is the following: While the Holder 
spaces are needed for a clean PDE theory, we want to prove a theorem whose 
natural hypothesis is about the L°° norm, rather than having to introduce an 



unnatural smallness condition in a Holder norm (see the proof of Thm. \ 7.8\ ). 
In particular, when estimating the nonlinearity, we need to use a priori small- 
ness in the L°° norm to avoid the singularity at u = 0, but may need to work 
in a subset of C° that does not impose smallness of the norm. 

We will also use a variant of Lemma 15.41 that deals with certain weight func- 
tions: 

Lemma 5.8 (Weighted product estimate). Let f G C^(]—a,a[ — )• R) with 
/(O) = and let < h ^ C^ill^ — i- R) be a positive weight function with 
|Vln/i| bounded. Then 

IK/ o w)w/i||c"(R") < C\\w\\L^(Rn)\\wh\\cc^{n") 

where the constant C depends on the norm of f and \^\ah\. The same es- 
timate holds with space-time norms overTU^, provided \dtlnh\ is also bounded. 

The crucial issue here is that there is no second term in which the Holder 
norm would apply to the unweighted function. Specifically, we are after the 
following 

Corollary 5.9. Let h{x) := (coshs)"'' on the cigar manifold Ai. Under the 
same hypotheses as in Lemma \5. <Sl we have 

IK/ o u;)w/i||c"(A4t) <C\\w\\loo^Mt)\\'^Hc'-{Mt) ■ 

Proof of Lemma \5. <g| and Cor. \5.9\ . Clearly || (/ o w)wh\\Lao is estimated in the 
claimed way; so we only need to estimate the Holder quotient. We estimate, 
abbreviating w{xi) =: Wi and h{xi) =: hi, 

{Wi - W2){f{Wi)Wihi - f{w2)w2h2) = 

{f{wi)wi - f {w2)W2){Wihi - W2/I2) + {h2 - hi)wiW2{f{wi) - f{w2)) . 
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With no loss of generality we assume /12 < hi. At least if Wi 7^ W2, we can 
obtain 

f{Wi)Wihi - /(W2)W2^2 _ 

\xi — 

f{Wi)Wi - f{w2)w2 Wihi - W2h2 



Wi - W2 



\Xi - X2\ 



+ hiWiW2 



f{wi) - f{w2) exp[ln /i2 - In hi 



Wi - W2 



\Xi - X2\ 



We then get the estimate 
\f{wi)wihi - f{w2)w2h2 



< 



\Xi - X2\ 



sup 

, l''l<ll'ii'IU°= 



dr 



+ ||w/i|koo llw I 



\wihi - W2h2 
\Xi — X2\°' 

sup |/'(r)| 

l'"l<li«'llLoo 



min{ 1 , I In /i2 — In /ii I } 



\Xi — X2I 



< C||w||2,oo + Ca||w||ioo ||ti;/i||j;^oo min|(i °',\\'V\nh\\Lood^ "| . 

The estimate persists by continuity if wi = W2- If space-time norms are desired, 
the same estimate applies to time Holder quotients. This proves the lemma. 

For the corollary on the cigar manifold, we only have to replace \xi — X2I with 
d{xi,X2) > \si — S2I, according to (1281) . □ 



6. SCHAUDER ESTIMATES FOR THE HEAT EQUATION 



The results and methods in this chapter are basically well-known. However, 
we need to elaborate on some details since we are also relying on some less 
familiar weighted norms. 

Also the impact of low regularity inhomogeneities needs to be taken into ac- 
count cleanly to cope with applications to the quasilinear case. We will use 
these estimates in the next section to get local existence for parabolic equations 
on uniform manifolds. Again, this will be done along traditional lines; however 
to control the noncompactness of the domain by means of uniform estimates 
requires recalling the proof details from the traditional case for inspection and 
reference. 

Some comments concerning the semigroup point of view and the issue of big 
vs. little Holder spaces will be given in Rmk. 17.21 below; suffice it to say here 
that these issues are not of any real significance for our purposes. 

So we begin by studying the heat equation 

(41) Vt-Av = dlr{t,x) + dih\t,x) + c{t,x), v{Q,x) = va{x) 
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l'^llc"{R^) 



5 



and use the following lemma, which is basically standard; but need to imple- 
ment some modifications to the usual versions. Namely, in order to prove a 
regularization estimate, we will also work with space-time norms that deteri- 
orate as t — )■ 0: We define the space C"(Ry) to consist of those functions for 
which the norm 

ll'"llc"»(R") •= Ikllc^CR") := max< sup t"''^||'u||c"([t,T]xR"), ||M||Loo(Rn) I 

^ '-0<i<T J 

is finite. We prefer the first notation for typographical clarity even when 
referring to functions in the larger space Cf (RJ). These norms refiect the 
short-time scaling [[(^^^(j^n) < Ct~"/^||fo||L°o(R"), which is satisfied by the 
homogeneous heat equation. The L°° norm is included in the definition in 
order to salvage the algebra properties Utttill^tii-Rn) < C*||^|Ic'»(r")II'^IIc"(R") 
and also the stronger 

(42) Ikf llc-(Rj) < C (||m||c-{r^)II^IU- + Ikllio 

which is a straightforward consequence of f l55]) . The first component alone 
would not control the L°° norm as the function u{t,x) = Int shows. 

Lemma 6.1 (Holder smoothing of rough sources by the heat fiow). There 
exists a unique solution to the heat equation ( HTj) in the space C°'(R^). It 
satisfies 

||'y||c"(RJ) < C ( ||fo||c"(R") + ll/llc^CR^) 

-I- Ta ||6||ca(j^n) + T 2"||c||car(Rn)^ 

and 

ll'i^llc^CRJ) < C ( ||wo||c«(R") + ||/||c"(R5;) 

(44) i-c, 1 \ 

+ T~\\b\\Loc(R,n) + T^"2"||c||^oo(Rn)j , 

and the regularization estimate 

, , ll^llc'-fR") < ( ll^o|U°°{R") + ll/llc-fR") 

(45) ' ^ 1 , 

where the constant C may depend on an upper bound for T. 
Going through the proofs, we will also see: 

Corollary 6.2. The very same estimates as in Lemma 16. il hold in the half 
space with Dirichlet (or Neumann) boundary conditions. 

Proof of Lemma \6.1\ By superposition, the estimates are assembled from four 
cases, in each of which exactly one of the quantities /, 6, c, Vq is non-zero. The 
estimate (H3|) for vq is classical, see for instance [Mj, IV, (2.2), with the critical 
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part being (2.13) there. Likewise we can refer to ([21], IV, (2.1), and (2. 8), (2. 9) 
for the estimating Wt' — Aw^^ = f^^ , and obtain our case with v = df^w^K 

We also note that the arguments for U and c in (H3|) can be given as (simpler) 
modifications of the estimates for the thus concluding the proof of (H3|l . 
However, as we will use the strengthened version f H4|) . we will give these es- 
timates explicitly in a moment. At the same time, we will later refer to (l43ll 
in order to derive (H5l) from it. For the convenience of the reader, we will 
also redo the detailed calculations pertaining to (H3l)-(H5l) in a self-contained 
manner in an appendix, with no-pretense to novelty. 

By r, we denote the heat kernel: r(t, x) := (47rt)~"/^ exp[— |xp/4t]. 

We estimate the contributions of ||&||l°° to the alternate bound (jH]): 

In terms of ||6||£,oo, we can estimate the spatial Holder quotient as 



I I J Q o Ft^^ I I 

With z := {x — y)/\/AT and h{z) := —2z exp(— |2;p) = Vr(i, z), this becomes 
where ( = {x' — x)/ -\/4r. 

Now for I s I < 1 , we can estimate the inner integrand by the integrable quantity 
supBi(z) \Dh\; for |s| > 1, we can estimate the inner integral by 2 J \h{z)\ dz, 
and combining the two cases, we get [v]x,a < CT^^~°''^^'^\\b\\Loci^-^ny The esti- 
mate of II w II LOO is similar, but simpler, providing a T^/^ factor. 

For the time Holder quotient, we estimate 

\v{t',x)-v{t,x)\< [ [ \VT{t' -r,y)\\b{T,x-y)\dydT 
Jt Jr." 



+ [ [ \Vr{t' -T,y)-Vr{t-T,y)\\b{T,x~y)\dydT 
Jo Jiv^ 

pt p pt' — T 

<CiVf-Vi)\\b\\Loo + \\b\\Lo^ / / VTsis,y)ds dydr. 

Jo ^R" Jt-T 



We only have to estimate the last integral yet. Carrying out the y integration 
first gives 0{s~^^'^), then the s integration yields 



C / [(f - r)-i/2 - (t - r)-i/2] dr = C{t' - tf'^ 



Hence [v]t;a/2 < CT^^'^)l^b\\Loo . 

We estimate the contributions of ||c||loo to the alternate bound fH4l) : 
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Finally, the c-term in (jH]), can be estimated in a similar way, only with T 
instead of VF. For the spatial Holder quotient we obtain 0(T^^2'). The L°° 
estimate gives a factor T in a straightforward way. 

For the time Holder quotient we can easily estimate 

\v{t',x) - v{t,x)\/\\c\\Loo < C [ [ln(t' - r) - ln(t - r)] dr , 

Jo 

which gives already the desired result, except for an extra logarithmic term; 
and this result would be sufficient for our purposes. But for good measure, 
let's prove the claimed, optimal, estimate: 

\v{t',x) - v{t,x)\ < / / T{t' - T,y)\c{T,x - y)\dydT 
Jt Jr." 

+ / / \^it' -'r,y) -T{t-T,y)\\c{T,x-y)\dydT 

Jo JR" 

< {t' - t + I)\\c\\lo^ 

with 

(46) / / \^{t' -r,y)-T{t-r,y)\dydT . 

Jo Jn." 

Now for ti > to, we have 

rf+ \>rf+ \ I |2> 2 2n(lnti - Into) , , 

r(ti,l/) < F(to,2/) \y\ <^:= 7 7 ^0^1- 

ii — to 

Hence 

|r(ti,2/) - T{to,y)\dy 



n/2 



IS"" ^1 / .nlnti— Into \ Inti — Into 



where 7(m) := /q" e~*s"/^~^ ds is the incomplete gamma function, whose second 
argument | we suppress, and of which we only need that it is increasing with 
finite limit at +oo, and with bounded derivative. With t\ := t' — t, to := t — r, 
ti — to = t' — t =: d, and {t — r) / d =: a a.s a. new integration variable, we get 
from g6]) 

I = CdjJ (^7(f(s + l)ln^) -7(fsln^))rfs . 

As sup |7'| is finite, the integrand is bounded by Cln^^^ = O(-) as s — )■ oo. 
Therefore I < Cd{l + ln(l + ^)) < Cd + Cd\n{l + ^). From this, it follows 

[v]t;a/2 < (CTl-"/' + CsUprfl-"/2ln(l + ^))||c|Uoo(R.) . 
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The sup is taken on when d = cqT with cq determined by a transcendental 
equation, and the value of the sup is therefore CT^~°/^, which proves the 
estimate for the time Holder quotient. 

We have thus concluded the proof of (jH]) and now turn to fH5l) . 
For the homogeneous heat equation with initial data Vq, we have first 

(47) \v{t,x)\ < \\vo\\loc . 
Then we have 

(48) \v{t,x')-v{t,x)\<C^-^^^\\vo\\Lo. 
because of the estimate 

v{t,x') - v{t,x) = / (T{t,x' -y) ~T{t,x -y))vo{y)dy 

= / / {x — x) ■ 'VT{t, X — y + s{x' — x)) vo{y) ds dy , 
Jr" Jo 

\v{t, x') — v{t, x)\ < \x' — x\ \\vq\\l^ / \'VT{t^z)\dz. 

JR" 



and (H7|) combined provide for t > t, 

^,^\v{t,x') -v{t,x)\ . {^(\x'-x\Y-^ r'^l'^ \ 

Similarly we can prove, for t' > t > t, that 

t' -t 

(49) \v{t',x)-v{t,x)\<C—^\\vo\\L^ 

and combine it with (H71) to estimate the time Holder quotient in the same 
manner. 

Taking the supremum over r G ]0,T], we conclude that HfH^a^pj^n-) < ||fo||L°°- 

Now let's look at the estimate of {dt — A)v = S^,/'-' with initial data 0. We 
decompose v = Vi + V2 where 

[dt - A)v, = X[o,./2] (t) ^1 |t=o = , 

{dt - A)V2 = X[T/2,T]{t) dfjp V2\t=0 = = V2\t=r/2 ■ 

For V2, the classical Schauder estimate ( H3|) applies on the time interval [r/2, T] 
and gives 

T"''^||'y2||c«([r,r]xR") < T"^^||f2||c<»{[r/2,T]xR") < C't"^^ || / || C"([r/2,T] X R") • 

For vi, we note that 

/.t/2 /. 

Mi'r^x)= / '9|,r(fr-s,a;-2/)f^(s,2/)cii/cis . 

JO JR" 
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Hence 

r.r/2 



(ir)"/^||fi||c"([r,r]xR-) < C'll^^idr, OIU-CR") < / (i^-s)"l/(s,-)||L-(R")C?s 

^0 

< C'||/||L°°([0,r/2]xR") • 

In conclusion 



t"/^||v||c"([t,T]xR") < C'(^T""^^||/||c"([t/2,T]xR") + ||/||l°°(RJ) 



If we can also estimate ||w||j;,oo in like manner, taking the supremum over r 
establishes the / part of (H5l) immediately. And indeed, we can estimate 

v{t,x)= [ [ d^^T{t-s,x-y)P{s,y)dyds 

Jo JR" 

5jr(t -s,x- y){r{s, y) - r{s, x)) dy ds , 



'0 ~'R" 

\v{t,x)\ < f j^^ _ ^y/^^, e-\^-'^'/^^''^\s-'^/^\x - ylll/ll^. dyds 



<c ts-'^/'it-sr/'-'ds\\fra.<c\mc. 

Jo 



This concludes the proof of the / part of fl45l) . The impact of b and c is 
obtained in a similar manner. □ 

Proof of Cor. \6.S[ The proof for Dirichlet boundary conditions is essentially 
the same: In (11141) and similar estimates, we merely have to replace T{t — 
T,x — y) with r{t — T,x — y) — r{t — t,x — y*), where y* is the reflection of 
y at the hyperplane bounding the half space; and of course we integrate only 
over the half space instead of R". This half-space heat kernel obeys the same 
Gaussian estimates, and no further changes are needed. In fll2ip and similar 
equations, where we had conveniently put the y into F and the x — y into the 
right hand side term, we can again swap them and now get expressions under 
the integral like (^diT{t — T,x — y)—diT{t — T,x — y*)) {b^{T, y) — }f{T,y + x' — x)), 
supporting the same estimates. The integrand within fH6l) now becomes 

\T{t' -T,x-y)- r{t' -T,x-y*)- r{t - t,x - y) + r{t - t,x - y*)\ 

< |r(t' -r,x-y)- Tit -r,x-y)\ + |r(t' -T,x-y*)- T{t -r,x- y*)\ 

and we may enlarge the integration over all of R" and re-use the old estimates. 

The same applies to Neumann boundary conditions. □ 

7. Quantitative global well-posedness of the linear and 
nonlinear equations in holder spaces 



Consider a smooth function on a uniform manifold M. Its derivative is a 
section in the cotangent bundle and its second derivative is a section in a 
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suitable vector bundle. We consider differential equations 

(50) ut = F{t,x,u,Du,D'^u) 

governing functions u on M. The precise structure of F is not of importance 
in this section. We only need that there is such a formulation. In local coor- 
dinates, it reduces to an equation of the same type, but now with standard 
derivatives in R". 

We assume that the equation can be written in local coordinates (and with 
the Einstein summation convention being in effect) as 

(51) ut = d^^{r{x,u)) + d,{b\x,u)) + c{x,u) 

with nonlinearities defined for u in an open interval f/ C R. 

It is not difficult to check that this property is independent of our choice 
of local coordinates and the partition of unity. We say that the equation is 
uniformly parabolic if there exist positive constants A, A such that 

(52) X\i\'<durix,u)Uj<M^\' 

holds in such local coordinates as in Def. 15.11 for m in a compact subinterval 
V GG U, where A, A may depend on V. We call the equation uniformly 
analytic if in addition for any compact V G U, all k > and all multiindices /3 

(53) \d^^d^Pix,u)\ + \d^d^b\x,u)\ + |9^9^c(x,m)| < Cr'=+I'^l(fc + l/JD'^+l^l 

with r depending on V. Again this property is independent of the chosen 
points Xj and i? - up to changing C and r. 

To address fl5T|) . we will also be interested in linear equations, which can be 
written in local coordinates as 

(54) Ut - dlia'^u) - di{Uu) - cu = f + d^g' , u{x, 0) = Uoix) 
assuming that the coefficients a*-' are uniformly elliptic and continuous. 
Let us tackle the linear case first: 

Theorem 7.1 (Inhomogeneous linear parabolic flows on uniform manifolds). 

Suppose that uq G C"(M), < a < 1. Then there exists a unique weak 
solution u G C"([0,T] x M) to and it satisfies 

(55) ||m||c«(Mt) < C (IImoIIc^cm) + ||/||l°=(Mt) + ||5'IU°=(Mt)) • 
where the constant C depends on T,X,A,n, ||a*''||(7a, ||6*||i^oo, ||c||Loo(]y[y-). 

Likewise, a unique weak solution to flMl) exists in C"(Mj') for uq G L°°(M), 
and it satisfies the similar estimate (!6T!) below. 

In the homogeneous case (f = 0, g = 0), the equation satisfies a comparison 
principle (u > if Uq > 0) and, provided the a'^^,b\c are independent of 
time, defines an analytic semigroup on C"(M) in a sense made precise in the 
following Remark \T7^ The equation also defines an analytic (in the sense of 
the following remark) semigroup on C;,(M) (the space of bounded continuous 
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functions with the L°°(M) norm) as well, and it satisfies, for < t < T , the 
regularization estimate ||M(t)||c«(M) < C'^ "''^||^o||l°°(m)j where C may depend 
on T. 

Remark 7.2 (Analytic Semigroups), (a) By an analytic semigroup in 
we mean a family t G [0,oo[ — > S{t) of linear operators on satisfying 
S{ti + ^2) = S{ti)S{t2) for all ti,t2 > 0, which extends to a complex sector 
and is holomorphic in the interior of that sector, satisfies a bound \\S{t)\\ < 
(j^CRet^ and for which {t,x) i— )■ S{t)uo{x) is in C° with respect to t and x, 
including at t = 0. (As t 0, this implies S{t)uQ uq uniformly, but not 
necessarily in the norm as would be required for the standard definition of 
analytic semigroup, which includes strong continuity at t = 0.) 

The analogous definition applies for an analytic semigroup in Ch; the continuity 
requirement in this case is that the mapping (t,x) 1— )■ S(t)uo{x) is in Cb with 
respect to both variables, including t = 0. This still implies S{t)uo — 
pointwise (but not necessarily uniformly) as t 0+. 

(b) Indeed, strong continuity at t = cannot hold in the space because, 
for Uq not in the closure of C°° within , smoother functions like S{t)uo 
cannot converge to uq in the C^-norm. Our estimates do imply that \\S{t)uo — 
^o|U°° — )■ as t — 7- 0. Moreover, if uq G with (3 > a, then \\S{t)uQ — 
^o||c« —'■0 as t ^ 0, because by two independent easy estimates. 



\u{t,x')-uo{x')-u{t,x)+UQ{x)\ ^ . \ {\u{t)]^,p + [uo]^,p)\x' - X 

< mm ' 



l3-a 



[ 2||M(t) -Mo|Uoc/|a;'-x|" 
and this minimum tends to uniformly in space as t — )■ 0. 

(c) Alternatively we could work in the 'little Holder spaces' oC"', the closure of 
C°° in instead of C7" itself. The estimates are the same, and the changes 
in the proofs are marginal. In oC", the strong continuity of the semigroup 
is restored, because uq G oC" can be approximated by (j) & for arbitrary 
(3 G ]a, 1[, and using the standard estimate 

\\S{t)uo - Mollc" < ||5(t)0 - 0||c- + \\S{t){u^ - + ||mo - 0||c" . 

(d) These observations ensure that solutions are continuous up to t = and 
the initial condition is understood in the obvious sense, albeit (for the classical 
Holder spaces) not in the functional analytic sense of norm continuity. Ba- 
sically, these subtleties at t = are of no concern to us as we are interested 
in asymptotics for t —> 00; on a technical level, this is borne out by taking 
iterative discrete time steps, and the finite time flow map has all the continu- 
ity, and even smoothness, we want, as elaborated upon at the end of the proof 



of Thm. 7.1 



Proof of Thm. \ 7.1\ We turn to Equation (15^ on the uniform manifold, i.e., 

Ut - dfAa'^u) - di{Uu) - cu = f + dig' , u{x, 0) = Uo{x) , 
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where the data a'^^ ,b\c, f,Uo are Holder continuous, and introduce a small 
parameter 7 to be chosen later. By the Holder continuity of the inital data uq 
there exists 6 such that \a^^{x) — a^-'{y)\ < 7 if d{x,y) < 6. We may decrease 
6 such that 36 is still a radius of uniformity. 

We use a partition of unity vf subordinate to a locally finite covering with 
balls of diameter 6, with local coordinate maps xi from these balls into R". 
As the a^^ are almost constant on these balls, we can, after an affine change of 
coordinates, achieve that a*-' — 6^^ is L°°-small on these balls. This introduces 
constants that can be bounded in terms of the ellipticity parameters A, A. We 
decompose u by 

vi := (fjiu) o « = E/ mi^i ° Xi) ■ 

In order to avoid clogging up the calculation with coordinate maps we define 
Vi •= Vi ° XT^ '■— ° puUbacks of the partition of unity 

and of the function u to the respective coordinate patches in R". This way, in 
the equations below, the coefficient functions as well as the quantities indexed 
with / live on R^ (or, in the case of ui, a subset thereof), whereas the function 
u lives on M^. 

Then vi satisfies the initial condition vi{0) = rjiUo^i and 

dtvi - Avi = TjidtUi - ^{riiui) 

(56) = dl [{d^ - 5'')r]iui] + [Ur^iui - 2a'^ {d,r]i)ui + g\] 

+ cr]iui + fr]i - g\d,r]i) + [a'^d^iji - b\dir]i)]ui . 

On the other hand, given u, we can use this equation to define the vi, and 
hence a mapping T : u ^ YliVi'^^i) ° Xh for which we want to find a fixed 
point by means of the Banach fixed point theorem (i.e., contraction mapping 
principle) and Lemma [6.11 On the support of r]i, a*-' — 5'-' is small in the L°° 
norm, and otherwise all coefficients are bounded in C". The rji are bounded 
in the L°° norm, but their norm will be large for small S. Given u, (H3|) 
yields, for each /, 

||^z||c"(Rj) < C (jl'qiUo^iWc'^ +7ll^«^dlc« + - 5*^||c« ||'7z^i/||L°° 

(57) +Ty^\\r]i\\c2,4ia^^,b\c)\\c4ui\\c^ 

+ T^~t||?7,||c"||/||c- +T^\\Vi\\c^'"\\9'\\c'-^ , 

where we have used the fine algebra estimate fl38|) on the term with the highest 
derivatives, but the simpler estimate HmwUc-c < C||m||c7q II'?^I|c<» suffices for the 
other terms since their coefficients can be made small by choosing T small. 
Since ui is the restriction of vl'^^k ^^'^ coordinate patch, and at most 

a fixed number N of terms (given by the ball packing Lemma 15. 2p contributes 
to this sum for each I, we can estimate ||Mi||c" < A^C'<5||m||c"{Mt)- 
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We want this estimate fl57j) for later reference, but we can strengthen it imme- 
diately: namely, a similar estimate based on is possible with only the L°° 
norms of b, c instead of the Holder norm; and also only the L°° norms of f,g, 
and correspondingly smaller exponents for T. 

Similarly, given two sources u and u for fl56|) . with the same initial trace Mq, we 
get for the difference vi — vi of the solutions the same estimate with u replaced 
by M — n on the right hand side and Uq = 0, f = = g. Thus 

(58) 

< (1 + T"/^||?7i||c-(Rn))||t^i - fi||c«(R?,) 

and similar with products rjirik of cutoff functions. Recall that ||Tm||c"(Mt) • = 
suPfc II l^z ^z'i^;%||c°(R5,)- Here, for each k, the sum consists of at most terms, 
where N describes the uniform local finiteness of the partition {fjf} as given 
by Lemma 15. 2[ 

Therefore, in order to show that the map u i— )■ Tu := '^{rjiVi) o;)(;; is a contrac- 
tion, with contraction constant < 1, we need to prove a similar contraction 
estimate (with smaller contraction constant i)/N) for the maps u ^ vi. Choos- 
ing 7 = 1/{2CN) we determine the necessary 6 and obtain a (large but fixed) 
bound for sup; All terms in flFTl) that contain potentially large norms 

of T]i are multiplied by a power of T, which we choose small to compensate 
for the norm of the cutoff {uq doesn't occur in the difference, and we again 
note that \\u — m||loo(Mt) ^ T°/^||-u — m||c«(Mt) because {u — u)\t=o = 0). So 
we conclude, for T sufficiently small (independent of Uq), that T is a contrac- 
tion (globally on all of C"(Mr)). An iterative procedure gives existence and 
uniqueness. 

An estimate 

(59) ||m||c"(Mt) < C (||mo||c"(M) + ||/||c"(Mt) + Ikllc^CMr)) 

similar to ( 155|) in Theorem 17.11 holds because Equation ( [57|) implies 

lirnllc.(M^) < iMc^^Mr) + CNT^'^CsWfWc^^ + CNT'2Cs\\g\\c^ + 

(oU) 

+ CN{Cs\\uo\\c^ + Cs\\u\\l^ + T^/'^Cs\\u\\cc.) , 

with constants now depending on norms of a*-^, 6*, c. We again estimate ||m||loo < 
IkolU- + T"/2||^^ _ uqWcc. < (1 + T"/2)||^iQ||^, + T"/2||w||^, and solve for 
IITwIlc-a = IImIIcc Likewise, by using the strengthened version of ( 157|) . based 
on dill), we eet (1551). 



Next, we prove the regularization estimate. Using fH^ instead of f H5]) . we 
can construct a local solution by Banach's fixed point theorem as before, with 
the analog of (1571) involving ||?7zMo||loo(m) < ||'"o||l°°{m) instead of ||?7ztto||c"(M)5 
and otherwise || ■ ||^a(]v[j,) instead of || ■ ||c"(Mt)- A slight modification of (158|) 
is needed: we use that ||'7z|Ic;<i(r") ^ maxjl, T"/^||?7z||(7a(j^n)}. Similarly, in 
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the analog of fl571) . the term \\a^^ — 5*-'||^q(r") II IU°° be estimated by 
max{7,T"/2||aU - 5*^'||c7«(Rn)} ||r/iM,||^<.(R„). 

This argument proves the short term existence and uniqueness for L°° initial 
data; it also provides, as in fIBU]) . an analog to flS^ : 

(61) ||'"|Ic'"{Mt) ^ ^ (||'"o||l°°{M) + ||/|Ic"(Mt) + llfi'llc"^(MT)) ' 

from which ||'u(t)||c«(M) < C'^^ "''^||wo||l°o(m) follows immediately when / = 0, 
^? = 0. 

Let us briefly note that these same estimates can also be used for differences 
Un ~ Um of solutions with, say, the same initial data uq and with f = = g, 
but different coeeficients a^ , 6^, c„ and a^, 6^, Cm respectively. They will then 
provide in analogy to (1571) that 

\\vin — Vim\\c"(Rr^) < 

< C (^-f\\r]i{uin - Uim)\\c'- + IK^ - (5*-'||c'^||77i(Mi„ - Ulm)\\L^ 
+ T^/'^\\r]i\\c2.^\\{a'^ ,bl^, Cn)\\c''\\Uin - Uimllc" 

+ 2j\\r]iuim\\c" + IK - all\\c^\\r]iUim\\L°- 

+ ^^^^II'7«I|C2.«||«^ - ali,K - bln,Cn - Cm)||c«||M«m||c") • 

Then, our estimate, applied to the respective fixed points m„ and Um, provides 
continuous dependence of the solution on the coefficients: 

(62) \\Un - Mm||c«(R5,) - ^'11^™ II C°(RJ) ||«^ " ClmyK " C„ - Cm)||c" • 

The comparison statement (nonnegative initial data give nonnegative solu- 
tions) would follow by commuting multiplication and differentiation and clas- 
sical maximum principles if the coefficients were twice continuously differen- 
tiable, and if the manifold were compact. The first obstacle can easily be 
dealt with by regularization: Approximate the a*-', 6* uniformly by smooth 
' K "with bounded norms, and such that the a^ , 6^ converge to a*-' , 6* in 
a (weaker) norm, with < P < a. The corresponding solutions Un then 
converge to u in the norm by the continuous dependence estimate (1521) . 

So we now write, for smooth coefficients in local coordinates 

Ut — a^^d'^jU — UdiU — cu = . 

Replacing u by e~^^u if necessary, we may assume that c < 0. Let p be the 
radius function constructed in Lemma [5.3[ Then —e{fit + p) is a nonpositive 
subsolution, provided fi is sufficiently big. Hence, if Mq is nonnegative and 
f = 0, g = 0, then u + e{fit + p) is a supersolution which tends to cxd as a; — )■ oo 
and which is nonnegative at t = 0. By the maximum principle it does not 
assume its infimum for < t < T and hence it is nonnegative. We let e — ?■ 
to arrive at the desired conclusion. 
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The construction for short time T can be iterated to get global existence in 
time, and hence a semigroup. The fact that we have an analytic semigroup 
follows from a standard argument: We introduce a complex parameter r and 
study, in a space of complex functions. 



The previous arguments work for all r > 0. We rewrite the equation as 



ut - ro(9j(a^^M) - ^^{b'u) - cu) = (r - tq) (9j(a^%) - di{b'u) - cu) . 



The right hand side is analytic in r. We obtain an analytic dependence on 
r, which is equivalent to having an analytic semigroup - as here we don't 
require norm continuity at time t = 0. Specifically, the same estimates we 
have made for r = 1 apply to the present case, for r G C, as long as \t — tq\/tq 
is sufficiently small (with the right hand side not going into /, but distributing 
over all the other terms in (1571) ). So our estimates extend to a small sector in 
the complex time-plane. 

Our estimates using the space for initial data will swiftly imply that 
the solution of (1541) defines an analytic semigroup on Cb(M), again in the 
sense of Remark 17.21 (We make the statement for Cf, rather than L°° to avoid 
technicalities in which sense initial data are taken on, technicalities that would 
not pertain to the focus of our problem.) 

To see this, let us first assume that our initial data Uq are uniformly continuous. 
Then they can be approximated uniformly by C" initial data wjf^. Using the 
estimate fl6T|) . we infer that the corresponding solutions u^'^^ converge in and 
hence uniformly to a (the) solution u for initial dat Uq. But the m^*^' themselves 
were continuous (actually they are in the unmodified C"(WIt)), so their limit 
is (uniformly) continuous. 

Non-uniform continuity of Uq can be handled by noting that for any continu- 
ous weight function w vanishing at infinity, wuq will be uniformly continuous 
if Uq G CbCM.) is continuous. Such a weight function can be constructed in a 
smooth manner based on Lemma 15.31 and then rewriting (jSlD for the corre- 
sponding wuq results in an equation of the same type, so that the result with 
uniformly continuous data for the conjugated equation gives the continuity 
result for the original equation with (not necessarily uniformly) continuous 
data Uq. The boundedness is of course maintained from the argument with 
the unconjugated equation. (In the specific situation M = that interests 
us for the fast diffusion equation, we can take w = (coshs)^''.) □ 

Basically the same proof yields 

Corollary 7.3 (Variant for Dirichlet boundary data). For equation (1541) on 
a smoothly bounded domain on M, with Dirichlet boundary conditions u = 
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(63) 





and initial data, there exists a unique solution in C", and estimate fl55|) 
continues to apply, as well as all the other conclusions from Thm. \ 7.1 



Now we can tackle the short-term estimates for the nonhnear equation by 
means of a shght modification of the proof of the hnear theorem: 

Theorem 7.4 (Short time smooth dependence on data). Suppose equation 
( l50l) can be written in local coordinates as (15T1) . namely 

ut = dl{r{x,u)) + d,{h\x,u)) + c{x,u) 

with nonlinearities defined for u in an open interval f/ C R, and that this 
equation is uniformly both parabolic and analytic. Let < a < 1, and choose 
two subintervals V CC W CC U . There exists T > depending on V and W , 
||^io||c'^{M); and the nonlinearities of the operator, such that for uq G C"(M) 
with Mo(M) C V there exists a unique weak solution u G C"(Mr) to flSTj) . 
Its values are in W . For every uq G C°(M) with range in V , there exists an 
R-ball in C"(M) and a T such that the map 

C"(M) D Br{uo) Buo^ue C"(Mr) 

is analytic. 

As stated, as of the present theorem, the existence time might depend on the 
Holder norm of the initial data. However, we will see soon that this is not the 
case and that the existence time is influenced by uq only through its supremum 
norm. 

Proof. The first step for the nonlinear problem, namely Equation ( 15T|) . 

= d^^{f'^{x,u)) + d,{b'{x,u)) + c{x,u) 

requires few changes. We rewrite the equation again as fixed point problem, 
this time in an ball about uq, within C". We choose Eq such that a 2eo 
neighbourhood of V still lies in W. We search a small solution in the form 
u = Uq + u, where u satisfies an equation of the same type; llMHioo will be 
assumed to be < 2^0 • 

Just as in the linear case, the equation for u can be localized using the partition 
of unity with 

vi = (fjiu) o xi^ , M = E/ Vii-vi o Xi) ■ 
Again ui := uoxY^ and rji := fjiox^^- One easily computes in local coordinates 

(64) dm - Avi = ^Ji^ + dibi + ci 

(with initial data 0) where, by a slight abuse of notation. 

If = Vi uo + u) - S'^u) , 

bi = riib\x, uo + u)- 2{djr]i) f'^{x, uq + u) , 
ci = r]ic{x, Uo + u) - (diTji) b'{x, Uq + u) + {df^rji) f'^{x, Uq + u) . 
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Using the considerations for the heat equation fHTj) . we again aim to construct 
the solution up to short time T, by means of the Banach fixed point theorem 
apphed in a closed subset S := {u \ \\u\\l°° < 2eQ , WuWc^ < M} of C"(Mt). 
Here M is a possibly large constant that will be chosen shortly and will depend 
on ||mo||c°(m)- The mapping is again T : m H- J2iiVi'^i) ° X;"^- ^^^^ ^ 
target contraction constant G ]|, ![• We also commit a priori to an upper 
bound To for T and calculate ||T0||c7a(Mj,) < ||TO||c«(Mto) =• ^- Then the 
choice M := A/{1 — i)) will turn out to be expedient. 

The estimates for the heat equation from Lemma 16.11 imply: 

(65) ||t;Hlc-(Rj) < C (llfiWco. + T'/^fbiWcc^ + ri-°/lQ||c<.) 
and for differences of solutions 

(66) \\vi - viWc^ < C {\\~fi - fiWc^ + T^'^hi - hiWc^ + T^-^/^Q - h\\c^^ 
with 

(67) fi-fi = r]i[ [duPix, uo + ^ + s{u- ^)) - 6''] ds (m - ^) . 

Jo 

The constants C in fl65l) . fl66l) are uniform as long as Mq, Uq + Uq + u have 
range within W . All norms in fl65|) . fl66l) are C°(RJ) norms. 

Smallness of the oscillation of duf^-^ over all pertinent arguments (x, u) G 
M X is crucial for the contraction estimate. Given 7 > 0, there exist 5, £1 
such that 

i9„r(x,M)-9„r(y,n)i <7 

if d{x,y) < 6 and — -ul < ei. By an affine change of local coordinates in a 
ball, we can actually assume S^^ to be the value of duf"^ at some {x,u). This 
will make our constants dependent on the ellipticity parameters A, A, but does 
otherwise not affect the estimates. We choose 7 < 2^ with bounding the 
number of (5-balls that can be packed in a 35-ball, a quantity independent of 6 
by Lemma Em Then, (1671) implies, with K a bound for \\duf''mL°°{MxW) and 
II^m/*'' IU°°(MxVK), that 

Wfl - /«||C"(RJ) < iWVli^l -^t/)llc-(R5,) + 

(68) +maXse[o,i] \\duP{x,uo + stt + (1 - s)u) - S'^\\c'-\\i]i{ui - ui)\\l-° 
< (^7 + K\\uo + sn + (1 - s)ii\\c- + l)CT"/^^ \\u - ?i||c-(M) , 

provided the partition of unity is made with balls of diameter < 6. With 6 
thus fixed, we have bounds for the C° norms of rji and its derivatives and can 
use the smallness of T, dependent also on Ar(||Mo||c'"(M) + M), to ensure from 
( l66l) and ( l68i) that T : u ^ ° X^'^ is a contraction from S fl C°(Mt) 

to C^iNlx) with contraction constant ^. 
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From (165|1 . we get a fixed (possibly large) bound for Tu := || '^(iliVi) o^i ^||c") 
and again the smallness of T ensures that HTmUloo < 2eo- We also estimate 

WTuWc^imr) < \\Tu - rollc^(M^) + ||ro||c^(M^) <i^M + A = M . 

This guarantees the applicability of Banach's fixed point theorem, and hence 
warrants local existence. 

The same solution can also be obtained by means of the implicit function 
theorem, because Thm. 17.11 guarantees a bounded inverse for the linearization 
operator. The analytic dependence of the equation on the data then leads to 
the analytic dependence of the solution on the initial data and parameters. □ 



We will want to control the deviation of the nonlinear flow from the linear 
flow, at least for short times. Contingent upon second order differentiability 
of the flow in the appropriate function space C°, this deviation should be of 
quadratic order in the norm of the space. But we will need an estimate that 
involves the weaker L°°-norm. So we claim: 

Lemma 7.5 (Linear approximation of nonlinear semifiow). Let w solve the 
homogeneous linear equation fl5^ . namely {dt — h)w = for initial data Wq, 
where in local coordinates, Lw = dfj{a^^w) + diiUw) + cw. Let w solve the 
quasilinear equation {dt — h)w = L{f{w)w) for the same initial data Wq, where 
in local coordinates L = L + 9jo6* + c and f is a smooth function from 
an interval about into R satisfying /(O) = 0. Assume the coefficients are 
smooth. 

Then, for sufficiently short time T and sufficiently small ||wo||l°° (dependent 



on the same quantities as in Thm. 7.4), there exists a constant K (uniform as 



T — )• Oj such that we have the estimate 

(69) \\W - w||c"(Mt) < K\\w\\c^(M.^) \\w\\loo(m.t) 

and from it the time-step estimate (with a different K): 

(70) 11^(2") - w(T)||c-(M) < -^lkollc-(M) IkolU-'CM) • 



Proof. The proof follows the Banach fixed point argument used in proving 
Theorems 17. II and 17.41 We decompose, as before, the solution wasw = ^rjiwi 
with wi = rjiw and {rjf} a partition of unity. Similarly, we decompose w = 
^ rjiwi in the same way as w, and refer to the proof of Thm. 17. H specifically 
Eqn. dnH). 

Now from Equation fl5B]) . we can copy 

^^^^ dtwi - Awi = dl [{d' - 5'')7]m] + [Ur^m - 2d^{d,7^i)wi] 

+ C7]iwi + [a'^{dl7]i) - b'{dir]i)]wi 
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and likewise from f l6^ 

dtwi - Awi = dl [{a^^ - 6'')r]m] + [Vr^m - 2a'^ {djTii)wi] 
+ cr]iwi + [a'^{dlr]i) - h\dir]i)\wi 

(72) + dl [d^r^if{w{)wi] + [U7^if{w{)wi - 2a'^{d,7^{)f{wi)wi] 
+ C7]if{wi)wi + [a'^idl7]i) - b\di7]i)]f{wi)wi 

+ di[b'r]if{wi)wi] + [crji - b\dirii)]f{wi)wi . 

Applying the heat equation estimate to the difference (which has initial data 0), 
we conclude that 

\\m - Wi||c«(Rj) < c(^j\\w - w||c°(Mt) + 7ll^«l|c"||^i^ - wWloo 

+ ||a*-^' — (J^-'llc" \\w — w||loo + T^/^||?7j||c72,q \\w — wWc^^ 

+ C\\ (a^^ b\ c, b\ d)\\co^ \\vi\\c^,c, \\fiw)w\\c^ . 

(All norms of w on the right hand side are space-time norms, even if notation- 
ally suppressed for conciseness.) Combining contributions from all coordinate 
patches, and using the contraction estimate for small times, we conclude that 

\\w - wWc-iMT) < K\\f{w)w\\c^(^MT) < K\\w\\c<-^Mt) \\w\\l°°(Mt) , 

i.e. (p]). 

From the maximum principle, we can estimate, similarly as in the proof of 
Thm. 17.11 briefly postponing details, that ||w||loo(]v[t) ^ C'||wo||l°o(m)- Now 
we assume ||tyo||L°° < l/(2-ft'C), and we obtain, still in the space-time norms, 
ll""^!! < ll""^!! + \\w — w\\ < \\w\\ + l\\w\\, hence ||w||ca(MT) ^ 2||iZ;||ca(MT)- From 
Theorem 17. we know that ||w||c'»(Mt) ^ C'||'^o||c«(m)- Combining these 
estimates, we conclude, with a new constant K' = 2KC, 

\\w — ty||c"{MT) < -^'||'W^o||c"(M) ||'W^o||l°°{m) , 
and this implies (170|) immediately. 

Let us now clarify the details of the comparison principle argument: After 
commuting the smooth coefficients in front of the derivatives, w satisfies (with 
appropriate new b\ c, b\ c) 

(73) wt - a'^d^iw + f{w)w) - UdiW - bA{f{w)w) - (c + cf{w))w = . 

We assume w w + f{w)w =: g{w) to be monotonic on an interval \w\ < 25i 
and choose 5q < 5 < 5q. Without loss of generality we 

truncate / smoothly to a constant outside the interval [— 25i,25i], leaving 
/ unchanged on [— so that the nonlinearity is globally defined and 
the equation with the modified nonlinearity is still uniformly parabolic. The 
modification of the nonlinearity will not affect solutions whose range stays in 
[— A smallness assumption on T dependent on 5i/5q will be imposed. 
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Namely, with max[_2s^,2Si] \f{'w) \ ='■ M, we want A > to be an upper bound 
for |c| + M\c\. For consistency, we require AT < ln((5i/(5o). 

We now compare w with w := e'^*(— ||iyo||L°° — ^{fJ't + p{x)), where p is the 
smooth radial function from Lemma 15.31 and A > as given above, and fi, e 
are positive constants. Given e > 0, we are assured that w < w outside a 
compact set, and also for t = 0. We need to show that the same operator as 
in (I73p applied to w is < 0. Then Thm. 12 in Sec 3.7 of Protter- Weinberger 
PP] guarantees w > w. Indeed, 

Wt — a^^d^jg{w) — IfdiW — b''-dig{w) — (c + cf{w))w = 

= Wt — g' {w)djjW — g" {w)a''^ {diw){djw) — (6* + Ug\w))diW — (c + cf{w 

= e^* ([A - (c + cf{w))]{-\\w4 - e{pt + p{x))) ~ep± eO(l)) . 

By choosing /i larger than the constant in 0(1), independent of e, this quan- 
tity is indeed < 0. Now we can let e — > and conclude w > — e'^*||^^;o||L°°• 
A similar argument can be made with w := e'^*(||w7o|U°° + £{pt + p(x)), with 
the inequalities reversed and an upper bound proved. □ 

We now return to the fast diffusion equation formulated in the relative L°° 
norm on the cigar manifold Ai, which is the motivating example for the equa- 
tions studied so far in this chapter. 

In order to get a priori control over the behavior of the dynamics for long 
time in the nonlinear case, we use a comparison principle with Barenblatt 
solutions to gain such control in the norm, and a parabolic Nash-Moser- 
DeGiorgi result to upgrade this control to a Holder norm. The a priori control 
gained from these arguments will allow to iterate the short-term existence from 
Thm. 17.41 and gain global well-posedness. 

It is of course well known that the initial value problem for the fast diffusion 
equation has unique solutions, and hence the initial value problem for v is 
well-posed in suitable function spaces. We need however a more precise result 
that establishes, in particular, differentiable dependence of the semiflow on its 
initial data. We now prove this in the Holder spaces C", remarking afterwards 
that the same proof extends to smoother spaces C'^'" and to spaces with 
more restrictive weights t] < (but not with the more permissive weights 
?7 > 0, which will be discussed in Chapter [TT]) . 

Lemma 7.6 (Relative bounds; cf [SJ). Suppose = mg < m E ]0,1[ 
and the initial data Uq G C(R") of a solution u to @ satisfies c^^Ug{x) < 
uo{x) < cUb{x) for all x G R" and some constant c > 1. Then there exists 
C = C{n,m,c, B) G [1, C)o[ such that (t, x) G [0, oo[ x implies 

(74) C-^Usix) < u{t, x) < Cus{x). 
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Proof. This follows from the same comparison argument used in Vazquez 
Thm. 21.1]. Consider a quotient of Barenblatt solutions: 

,75) f^f4^y"-"V(i,2,,)-.B.(l + 2p.)- + |,P 



UB{y) J B+\yY 

Given any number c > 1 we can make the quotient Pb+{'^, y)/uB{y) larger than 
c (equivalently, the right hand side of (175|) sufficiently small) by first choosing 
T large and next 5+ > small. This yields the inequality 

pB+{r+,y) > cuBiv) > uo{y) 

for some 5+ > small enough and r+ large enough. By the comparison prin- 
ciple, the solution of ([T]) with initial condition po '■= Uq satisfies p(r, y) < 
Pb+{t + 'T+^y) < C+pB^{T,y) for all (t, t/) G [0, oo[ x and some exphc- 
itly computable C+ = C+{n,m,T^) < oo. The corresponding transformed 
solutions to equation (jl]) satisfy u(t,x) < C+Ub_^_{x). Since 1 < Ub^Iub < 
{B^/ B)^^^^^~^^ this completes the proof of the second assertion in fTMj) . 

The same estimate with large and r_ > — l/2p sufficiently negative can 
be used to get p(r, y) > p^ {t + T_,y) > C^Pb_ (r, y) and complete the proof 
that u{t,x) > Ub{x)/C for some C G [l,oo[ and all {t,x) E [0, oo[ x R*^. □ 

In order to use the DeGiorgi, Nash and Moser result, we note that in local 
coordinates we can write equation flSlj) in divergence form as 

ut - di{{duP)ix,u)dju) - [{duV){x,u) + {diduP){x,u))djU = c{x,u) 

where 

c{x,u) = c{x,u) + id,b'){x,u) + (9jr^Kx,n) 

is bounded if u is, i.e., if u(Mt) C W. Let us clarify the notation here: 
expressions of the form di{g{x,u)) refer to a partial with respect to Xi in all 
locations including implicit in u, whereas {dig){x,u) would refer only to the 
explicit occurrence of Xi in the first argument of g. 

We consider this equation as a linear equation of the type 
(76) Ut- di{a!'^dju) -UdiU = c 

where a*-', U and c are bounded measurable functions and there exist A, A > 
with 

in the sense of quadratic forms. Now we can apply 

Theorem 7.7 (DeGiorgi-Nash- Moser) . There exists < < 1 depending 
only on A and A so that any bounded weak solution u to (!76|) in the cylin- 
der [0,2i?2] X B2R lies in C"o([i?2^ 2/?^] x Br) and 

\W\\c°'0{[R-^,2R^]xBb) < C(A, A, R, ||6||l°°) (||M||L°°([0,2/?2]x_B2fl) + l|c||L°°([0,2R2]xB2i{)) • 

45 



Proof and statement can be found in Ladyzenskaja, Solonnikov and Ural'ceva 
[31], III§10. 

We can now state: 

Theorem 7.8 (Long-time smooth dependence in of the FDE on data). 

Let a G ]0, 1[, m G ]0, 1[ andm > itlq = Then there exists R such that for 
vq G C" with \\vo — 1||loo(R") < R there exists a unique solution v G C"(A^t) 
to fl2T|) for each T. The map 

has continuous derivatives of all orders. It is bounded in the sense that < 
v{t,x) < C for all {t,x). 

Proof. The local existence part is an immediate consequence of Theorem I7.4[ 

We are now in a position to see that the local existence time asserted in 
Thm. 17.41 is controlled only by the norm of the initial data, and not by an 
otherwise conceivable deterioration of the Holder norms. 

Comparison of u with K ±{mt + eo) (with m large enough to control c, db and 
d'^f) shows that there exists Tq independent of the Holder norm of the initial 
data so that u cannot leave W before that time. We again turn to the local 
uniformly parabolic equation which we write as 

Ut = di{a^^{x, u)diu) + 6*(x, u)diU + c(x, u). 

i 

By Theorem 17.71 and uniformity of the manifold there exists ao > and C so 
that ||'u||c"o([T,min{T,To}]xM) < C with C depending only on W (via the bounds 
of derivative of the nonlinearities, and the ellipticity of the operator) for all 
r > 0. 

Localizing again on sufficiently small spatial and temporal scale we obtain 
again local equations 

■Uf - A-u = i^f^ix, uo + u) - S''^{uo + u)^ + diU{x, u^ + u) + c{x, uo + u). 

Now 

9j.(r(x, uo + u)- 6''u) = dlif'\x, uo + u)- r(x, uo) - 6''u) + d^f'^x, uo). 
Multiplying by a cutoff function we obtain 

\\u\\c'-([to+R2^to+2R2]xBR{xo)) < c\\u\\L'^(^[toA,+2R^]xB2R{xo)) 

which gives the bound of Thm. 17.41 as long as u assumes values in W. 

Therefore the only way that a solution may cease to exist is by v tending to 
zero or infinity. This is impossible by the pointwise bounds of Theorem 17.61 
This proves the global semifiow property. □ 
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Remark 7.9 (Smoother spaces). We can carry out the same proof in C'^'" 
spaces, using the analog of fH5]) for these spaces, or by looking at systems of 
equations for the derivatives. 

We observe that without any change in the argument we obtain a similar global 
result in the function space for f] < 0. This space is defined by the norm 

(77) ||/||c.:=||(cosh.)-V||c^. 

Let us be more precise. We define v = (coshs)~''f, and then i) satisfies an 
equation similar to f l2T|) . for which the same proof carries over. In this pro- 
cess, ?7 < is needed to guarantee that still embeds into L°°. Otherwise 
quadratic and higher order terms accumulate positive powers of cosh s and 
cannot be controlled anymore. 

We will also obtain certain analogs of Thm. 17.41 and Lemma 17.51 in the case 
rj > in Chapter [TTJ They are of a mixed nature inasmuch as the hypotheses 
and statements still need to enforce the boundedness of w that is not enforced 
by the norm itself. 

8. The spectrum of the linearized equation 

Since the spectral calculations of |19] and [20] have been done in a different 
framework, some explanations are at hand to obtain a valid comparison. Here 
we have adopted the normalization conventions of the announcement [T9] . 
which differ from those of [20j. 

Our fast diffusion equation ([1]), (jl]) differs from the one in [20] by a factor 
1/m in front of the Laplacian and a factor in front of the rescaling term. 

This difference is explained by a correspondence = jz^^ ^Ho] = ^' 

C[2o] = Yi^B. T coincides between [20] and here, but = 1/(3. The 
linearization operator in [20] was 

H[2o]^ = -m<-^A^ + X ■ V$ 

(of which we should drop the m in front of the Laplacian and add a factor 
of in front of the rescaling term to account for our space variable, and 
which was defined there as the negative linearization operator). Let us import 
the spectral results from [20], but adapted to the notation conventions of the 
present paper, which coincide with those of \T9\: 

Theorem 8.1 (Spectral theory of the linearized operator [20]). The operator 

(78) H : $ ^^ <-2v ■ [u^Vm = m^-^A^ —x ■ 

\ J B V B i B l-m 

is essentially self-adjoint on the Hilbert space defined by the norm (J V^'p dxY^"^ 
(with constant functions \E' modded out). The restriction o/H to the eigenspaces 
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of the spherical Laplacian As with eigenvalue + n — 2) (where i 
0,1,2, .. .) has continuous spectrum for 



A < A 



cont 



-i{i+n-2)-{ 1+1 



[l-mf ' \2 ' ' 'I 1-m 
and finitely many eigenvalues 

A^fc = —^{£ + 2k) + Ak{k + £ + I - 1) = -(^ + 2k)p - n£ - 4/c(l - I - k) 
for integers k satisfying 

0<^<HT^-f + l-^]=Mf + l-^], 

with corresponding eigenf unctions ipik{f^)Yi^{x/\x\) where the Yi^j_ are spherical 
harmonics and 

k + i- l-p/2, -k 

£ + n/2 '"5^ 

which are polynomials of degree I + 2k in r. 

Moreover, the operator L studied here is still not exactly the same as the op- 
erator H just imported from [T9], [20]. This is due to a different linearization. 
In [19j, [20j, the linearization was defined in terms of the pushforward of a 
measure, (/ + £:V\l/)#Us, whereas here we have the linearization + ev). 
At least on a formal level, 

((/ + eV^)#Us) (x) = Us ((/ + £V^)-^a;) / det(J + eV^) = 

= {ub{x) -eVugi^x) ■ V^(a3))(l - e trace L'(V^)) + 0{e'^) 
= Ub{x) - e{VuB{x) ■ V^(a3) + Ub{x)A-^{x)) + 0{e^) 

= {ub-bV ■[uBV^]+0{e^)){x), 

so —V = Mg^V • (mbV\1') =: A\E'. So we should have L o A = A o H. This 
can be confirmed by a look at the 'factorized' forms in flTSl) and fl23|) . Since 
A = M^""^ o H, there is also an easier conjugacy: L o u^""* = u]^~^ o H. (In the 
notation fl25l) . this means H = L^=„2-) 

In Theorem 18. 11 an eigenvalue Aqo with eigenfunction a constant, does, strictly 
speaking, not exist, because constants were modded out in the Hilbert space 
^UB defined by the norm (J Us| V\l/p(ia;)^''^. In our context however, the con- 
stant function t/^qo (which becomes proportional to u\~"^ under conjugacy) will 
correspond to the derivative of the Barenblatt with respect to mass: indeed, 
from {dBUB)/uB = yz^^b"'"- The mass was of course fixed a priori in the 
(mass-transport based) linearization formalism of [19] [20]. 

We can also introduce the Hilbert space L^2-m, which is defined by the norm 

Ug 

(/ u^^l'^l'^dxY^'^ . Now first by abstract functional analysis, then by explicit 
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integration by parts, we have 

B B Ug 

Thus, as the quartet [8] also found, we get a Hilbert space isomorphism H^/^ : 
W^'^ — )■ L'^2-,n and H is self-adjoint with the same spectrum and eigenfunctions 

in L^2-m as well, except that Aqo = with the constant eigenfunction is genuine 

on L^a-m- 

Ug 

For L, which is self-adjoint in the Hilbert this means (at a formal 

level, since we want to consider L in a Holder space setting) that the eigen- 
functions for A^A; are now v}~^il)^k{''^)Y^^{x/\x\). They will lie in the Holder 
space only if they are bounded. Different weights in the Holder space can 
however be used to incorporate these formal eigenfunctions into the space, or, 
equivalently, we can conjugate the operators by an appropriate weight and an- 
alyze the conjugated operators in the unweighted space C". In any case, the 
critical growth threshold for a function to be in L^m differs from the growth 
threshold for a function to be in C". This affects the spectral theory. We 
define 

C'^^-{M) := {g := (coshs)V I / G C^^^{M)} 
^ ' ||5'llc^fc.-(-A^) := ll(coshs) ''5(||cfc.<.(Ai) . 

Also note that rjcr = \ — ^, combined with ( 1251) and (I5U]) . yields the isometry 

II/IIl^^(R") = ||(coshs)-''='-/l|L^(A4) =: II/IIl2^,(a^) 

and thus displays that rjcr gives the critical growth for selfadjointness on M.. 

Given an unbounded operator L : domL — > defined on a linear subspace 
domh C C", recall that A G C is said to be in the spectrum of L if (L — XI) 
does not have a bounded inverse on C". There are two subclassifications of 
spectra that we also refer to: into point, continuous, and residual on the one 
hand, and into essential and inessential on the other. More specifically, A is 
said to be point spectrum if L — A/ fails to be injective, continuous spectrum if 
L— A/ is injective and its range is dense but not closed (in which case (L— A/)~^ 
fails to be bounded), and residual spectrum if L — A/ is injective but its range 
fails to be dense. We say L — A/ is essential spectrum if range(L — XI) fails 
to be closed (finiteness of its codimension, when closed, and of the dimension 
of its kernel, always being satisfied in the examples below; c.f. [301 IV §5.6]). 

The calculations from [20] can be reused, to give the following result: 

Theorem 8.2 (Spectrum of L in Holder spaces). Given r/ G R, the formula 
Lv := Ug^'V ■ [ub'V{u^~^v)] defines an operator in the space C^{Ai) with do- 
main C^'°'{Ai). Let Li be its restriction to the eigenspace of As for eigenvalue 
-£{£ + n-2) (where i = 0,1,2,.. .), i.e., L(/(r)F,^(u;)) = (L,/(r))r,^(a;). 
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The domain ofLe is 



C, 



2,a 



considered as an operator in 



/(0) = /'(0) = /"(0)=0 ifl>i 

/(0) = /'(0)=0 j/^ = 2 

/(0) = /"(0)=0 if 1=1 

/'(0) = ifl = Q 



CI, := {/gC"[0,oo[|/(0) = Oz/^^0}. 



Let r]f,T = \ — 1, the critical value for rj introduced after fl26l) . The spectrum of 
he consists of finitely many eigenvalues, plus an unbounded connected compo- 
nent consisting of a filled-in parabola (which degenerates to a ray if rj = rjcr)- 
The boundary of this parabola forms the essential spectrum. More precisely: 

(1) If rj = rjcr, then he has only point spectrum. This consists of essential 
spectrum for A < A^"'^* ( as in Thm. \8.1\) in addition to discrete eigenvalues Xek 
of finite multiplicity for integers 0</c<|[| + l — £]. 

(2) If rj > rjcr, the spectrum is still only point, with the discrete eigenvalues 
exactly the Xek for < k < + 1 — ^— 1^] — ?7cr|] dnd the parabolic region 



ReA < 



ImA/2\2 
rj-rjcr^ 



i.e., ReA/T<?7 + l- | = |?7- rjcrl, where T = + n - 2) + (| + 1)^ + A . 
The largest ReA occurring in the essential spectrum is precisely at Xf^^ := 
iv — VcrY ~ (f + -'-)^ ~ ^(^ + n — 2). (This parabolic region includes now those 
eigenvalues Xik that, in comparison with the case rj = rjcr, have been 'lost' from 
the original list of eigenvalues.) For i = and rj = (unweighted C^), this 
threshold Aq°q* coincides with Aqi- 

(3) If rj < rjcr, the same formula Re Vt" < I?? — "ilcr\ characterizes now the 
residual spectrum (h — X has closed range with codimension 1 ), and the same 
formulas as in (2) apply for the eigenvalues Xgk, which now make up the only 
point spectrum. 

In cases (2) and (3), the boundary of the parabolic region, namely Re-y/r = 
?7 + l — | = |?7 — rjcr\, is essential spectrum (range L — A not closed). The 
remaining spectrum is non-essential: the kernel of h — X is finite dimensional, 
and its range is closed with finite codimension. 

The eigenf unction for eigenvalue Xek is vek{\x\)Ye^{x/\x\) where 



vek{r) = Ug ™V£fe(r) = (coshs) ?/'£fc(sinhs) 



/ sinh^ 



Vcosh 



X 2F1 



k + 



■p/2, -k 
i + n/2 



sinh^ , 
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T = £(£ + n-2) + (| + l)2 + A 



H in Hilbert space L a-m or W}, 
L in Hilbert space L^m 

® 



1,2 

B 



continuous 



parabolic region: Re a/t < \f] — f]c 



with T] > T], 




ponit 



with Tj = Tjc 



E _ 1 

2 



point 

with T) < Tjc 

residual 



® 



t ® 



Figure 2. The spectrum ofh^, in various spaces, schematically. 



Remark 8.3. The domain of L is dense in the little, but not the large. Holder 
spaces C". The latter fact seems a bit annoying from the point of view of 
abstract operator theory, but will not cause us any trouble. 



Remark 8.4 (Relevant parts of spectrum). We will rely only on the value of 
the essential spectral radius of the semigroup generated by (essential spectral 
abscissa of L^j, and the eigenvalues outside the essential spectral radius, but 
not on any other aspects of the spectrum. 

After decomposition into spherical harmonics, the value of the essential spectral 
radius could be obtained qualitatively by asymptotic methods (see, eg., Ch. 2 
of Fedoryuk [21]j without reference to explicit solutions in terms of special 
functions. The precise nature of the spectrum is instructive to know, but not 
relevant for our arguments. Note that in contrast to the weighted space used 
in [2U], in the present, unweighted space, the onset of the essential spectrum 
o/Lq is precisely at Aqi- 

Proof of Thm. \8.2[ Lest there be any doubt, we point out beforehand that 
the sign of t], rj^r will not be relevant for the proof of this theorem (but will 
become relevant later). The claim that L defines an operator in with the 
domain C^'" is equivalent to the claim that the conjugated operator = 
(coshs)"'' o L o (coshs)'' is a mapping C^'° — )■ C", and this is clear from the 
definition, see fl2S]) : or the equation in cartesian coordinates to see that there 
are actually no singularities at the origin. By slight abuse of notation, we let 
operate on functions /(s) rather than f {s)Yn^{x / \x\) . Recalling (!25|) . we 
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have 



L^,J7 :=(cosh s) oLi o (cosh s)^ 

= (ds + ^^""r^^ g. - (tanh e{i + n- 2) 
V sinh 2 s 

(81) / 2m \ , o 

+ [n h 2?7 ) tanh s Og 

\ 1 — m / 

/ \ ( 1 2m 4 \ ,2 

+ n(?7 + 2) + 7] tanh s . 

\ 1 — m 1 — m/ 

For the claimed domains C^'", notice that if /(r)Y£^(a;/|a;|) is in C^'"(A^), 
then / is C^'° by restriction to rays. Moreover, any C^'"(A^) function can 
be written, near the origin, in the form Cq + Ci ■ a; + c^r^ + h^ix) + o(r^), 
with /;,2 a harmonic polynomial of homogeneous degree 2. (Readers looking 
for details on harmonic polynomials can find them on pp. 159-163 of [5]). If 
this function is also of the form /(r)Y£^(a;/|a;|), then we can multiply with 
various Yf^' and integrate over the unit sphere. Unless £ = 0, choosing £' = 
gives = Co + C2r^ + o(r^), hence Cq = C2 = 0. Unless £ = 1, choosing = 1 
gives = Ci. Unless £ = 2, testing with = 2 gives /i2 = 0. The claims about 
/ and its derivatives at follow. 

Now we need to see, conversely, that when / is in the claimed domains, then 
fYiy^ is in C^'"(A4). It is only in a neighbourhood of the origin that this is 
nontrivial. Using Schauder theory, it suffices to show that A(/(r)Y£^(ci;) = 

( J^- + ^ ~ ^^^-7^--^)/('^)^/i(i^) is (uniformly) in a punctured neighbor- 
hood of the origin. Use of the initial conditions at in each case establishes 
this in a straightforward manner. 

The analysis of solutions to (L^ — A)?/^ = in terms of hypergeometric func- 
tions carried out in Sec. 4.2 of [20] carries over with only trivial modifications, 
and the discussion whether such solutions lie in C°(A^) reduces to growth 
at infinity (the smoothness being trivial). This establishes the precise point 
spectrum claimed for all cases, and the eigenfunctions Vi-^. 

We now assume we are not in the eigenvalue case. Translating the resolvent 
formula (4.36) from [20j, we can write any solution to (L^ — A)^ = w as 

(82) 

u(r) = v^ir) [Aq + wviUgV"' ^dr) +t>i(r) (Ai + J^wv2U^r^ ^dr) 
u'{r) = V2{r) (Aq + /J wviu^r^^^dr) + v[{r) i^Ai + j\^wv2u'^r'^~^dr) . 

As r — !■ oo, we have u'^r"'~^ ~ Q^-p/^-ricr Q^y^^^ 

V2{r) ~ cr''"'"^ , fi(r) ~ cr''"'+^ . 
As r — 0, we have u'^r^~^ ~ cr"^^^ and 

/ \ 2—n—f / \ f 

V2[r) ~ cr , Vi[r) ~ cr . 
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First assume Re Vt" > \v ^ VcA and let w = 0{r^) as r ^ oo. Then the 
second integrand converges near oo and we may specify its upper hmit as oo. 
A growth estimate for u requires now Ai = 0, because else the dominant term 
j\^^Vcr+VT jg |-]^g space. The dominant term as r — >■ is Agr^"""^, so we 

need = to get m(0) = in the case i > 1. In the case £ = 0, we still need 
Aq = for m'(0) to vanish. With these choices we do indeed get a solution u 
that satisfies the bounds near and infinity required for it to be in the domain 
of he (m"(0) = in the cases £ > 3 uses w{0) = 0); the smoothness is trivial. 
So these cases are not in the spectrum. 

Next we assume Re i/r < \t] — T]cr\- We are then in the case t] < rjcr because 
the other cases have been dealt with as point spectrum. If G range(L£ — 
A), then wviU^r"'~^ dr = by symmetry of with respect to the inner 
product in L'^(R~^ ,u'^r"'~^dr), which is defined in this particular case because 
the integrand is bounded by (9(^7-»?+R-eVT-»?cr-i^ ^^j. i^rge r. Let us assume w 
satisfies this orthogonality condition, which restricts w to a closed subspace 
of codimension 1. We also have the convergence that allows us to have the 
upper limit oo in the second integral of the resolvent formula f l5^ . Moreover, 
if w = (Le — X)u we need = in the resolvent formula for the same growth 
reasons. We can then write 

poo roo 

u{r) = AoV2{r) — V2{r) / wviu^r'^~^dr + Vi{r) / wv2U^r'^~^dr . 



We now need = 0, because the first term is the only one that would grow 
faster than O(r^). As before, we now confirm that u does have the required 
behavior as r — to be in the domain of — A. We have (non-essential) 
residual spectrum in this case (injective, range closed and with codimension 1). 

Let us now return to the point spectrum in Rca/t < \'ri — ricr\ in the case 
Tj > rjcr to see that this spectrum is non-essential. Indeed, we want to show 
that — A is onto in this case. We cannot choose oo as upper limit of 
integration in f l82p and choose 1 instead. Ai is arbitrary anyways, since Vi is 
a bona-fide eigenfunction. The choice = is required for the same reasons 
as before, and with this choice we do get a solution u in the domain for every 
choice of w in the space C^^. 

Finally, we consider the cases on the boundary of the parabolic region. Let 
first Re a/t = rj — rjcr > 0. We know we have point spectrum, and we want to 
show that the range is not closed. Not being assured of convergence at oo, we 
write 

-u(r) = V2{r) ^Aq + J wviUg'r'"'~^dr^ + Vi{r) (^Ai - wv2V^r''~^dr 

For w to be in the range, it is necessary that the function r ^ wv2U^r'^~^dr 
remains bounded as r — > oo, which is not automatic: w G would allow for 
logarithmic divergence. Assume this hypothesis is verified. Then the choice 

53 



Aq = does give a solution that satisfies the right bounds (and automatically 
the needed smoothness, too). With this characterization of the range, it can 
be seen that = r'^^'^~'^'=^^V2x{'i") / Inr (where x is a truncation at 0, constant 1 
for r > 1) is not in the range, but Wn = r'^^^~^''''^V2x{^') / i^^i" + j() is in the 
range, for every N, and Wn ^ w^, uniformly, and also in the Holder norm. 

Now, for T] < rjcr and Re a/t = rjcr — f], we can write the second integral as 
again, and need Ai = 0, Aq = 0. We get a similar conclusion, namely that w 
is in the range if and only if r h-t- wviu^r^~^dr is bounded as r — ?■ oo. With 
the same argument as before, the range is not closed. Is the range dense or 
not? We claim it is not dense, thus identifying the spectrum as residual rather 
than continuous. Indeed, if ||5'||c^ < then w(r) := r^^''~'''=''^'yi(r)x(r) + g{r) 
parametrizes an open set of functions that is not in the range, where the 
implied smallness of the weighted norm of g prevents it from compensating 
for the logarithmic divergence of the integral caused by the first term. □ 

The spectrum of L in is the same as the spectrum of (defined in fl25|) ) 
in C". So we can now study the differential equation vt — L^v = 0. It fol- 
lows from Rmk. 17.91 that the linear semigroup S^(t) generated is analytic on 
C°'{A4). We want to obtain the usual semigroup estimates on the comple- 
ment of the space spanned by the eigenfunctions vek- Major work goes into 
the boundedness result: || exp[tL^]|| < C exp[coo{ri)t] without an extra e. The 
loss of self-adjointness makes this result non-trivial, in particular since we have 
established non-zero Fredholm index for the onset of the essential spectrum 
of L^. We claim 

Theorem 8.5 (Semigroup Estimates). Equation f l2^ defines an analytic semi- 
group S^(t) = exptL^ on C"(A^) or on L°°(M), in the sense of Remark\7.2. 



The essential spectrum ofSn(t) is contained in 5(0, e'^°°('')*) and in no smaller 
ball, where c°° (77) is given by (1261) and equals the A™^* from Thm. \ 8.2[ Outside 
this ball there are only finitely many eigenvalues e^^* with Cooiji) < Ai < A2 . . . 
where the Xj are the Xek from Thm. \8.B 

The spectral projections onto the eigenspace corresponding to a set of eigen- 
values Xik of the operator L in L^m(R") = L^^^{A4) are, by the same formula, 
also well-defined mappings in C"{Ai), or correspondingly weighted L'^{A4), 
provided the L'^^^{M.)- eigenfunctions for Xek are still in the space C^{A4). 

Let Vj be the eigenfunctions corresponding to the Xj in L^{Ai). Then, for 
V 7^ Vcn each solution v = S^(t)-Uo with initial data Vq G L°° can be written in 
the form 

(83) v{t, a;) = ^ Cje^'^Vj{x) + Vres{t, x) 
with 

(84) sup e"''"^''^* Uteres < 00. 

t>i 
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For 1] = rjcr, we at least get the same conclusion with c°°{ri) + e instead of 
c°°{ri), for arbitrary e > 0. For rj = rjcr, we also get the analog of flMl) with 
the L'^{M.) norm (and no e). 



Proof. It has been proven in Theorem 17.11 that the semigroup is analytic. The 
spectral statements follow from Thm. I8.2l by the functional calculus for analytic 
semigroups (e.g., Cor. IV, 3.12 in Engel, Nagel [23]). In using this reference, 
note that strong continuity in is assumed there; so we get the estimate of 
the spectral radius for the little Holder spaces oC" directly, in view of our 
Remark 17.2( c). However, the result also carries over to the big Holder spaces 
C", because the semigroup maps initial data continuously into C" C oC" 
(for a' > a) in any short time to, as a consequence of the regularization 
estimate in Thm. 17.11 

Clearly the projections C°(A^) 3 u{s,u) i — > Ue^{s) = J u{s,u)Yifj^{u) du e 
C°[0, oo[ are continuous. If £ > 0, the functions ug^ satisfy U£^(0) = 0, and we 
denote the space of those C°[0, cxo [-functions that vanish at by Cq [0, oo[. 

Conversely, the imbeddings Co[0, oo[ 3 u{s) i — )■ n(s)Y£^(a;/|a;|) G C°(A^) 
are continuous for i > 0, and likewise the analogous imbedding C"[0, oo[ — ?■ 
C"(A^) in the case i = 0. The L°° estimate being trivial, we deal with the 
Holder quotient: 

\f{r,)YM-fir,)Yiu,)\ ^ \f{r,)-f{r,)\\YM\ ^ \f{r,)\\Y{u^)-Y{u2)\ 
d{{ri,uJi),{r2,UJ2))" - rf(ri,r2)" c?((ri, Wi), (ra, Wa))" 

The first term is bounded in terms of [/]„. The second term vanishes if £ = 0. 
If however £ > 1, we use that /(O) = and therefore \ f{r)\ < O{d{0,r)°') = 
0(s"). Then the second term is bounded by 

S^d{uJi,UJ2) 



max{|s2 — Si|", (tanhs2)" d{uji, W2)"} 

where Si > S2 without loss of generality. The estimate follows immediately. 
Finally, the inner product {v,V£k)L'^^ is well-defined for v G C"[0, oo[ provided 

7] + £ + 2k < p. By the constraint on k, we have indeed r] + i + 2k < i] + 
^ + I — \ri — rjcrl = p + rj — rjcr — If] — rjcrl < P- So again these projections are 
continuous in each C^, and so are trivially the corresponding imbeddings. 

So, since for given m < 1, there are only finitely many eigenvalues Xek, we 
conclude that the spectral projection Q onto the span of their eigenspaces that 
was naturally constructed within the L^m Hilbert space framework, namely 

{v;vekYe^)Li^ 
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defines, by the same formula, a spectral projection in the C^{M) framework, 
too. Here, k{e) = i[f + 1 - f - jr] - r/^rl] is from Thm. O Let P := 1 - Q be 
the complementary projection in C^. 

We turn to the semigroup estimate and consider a solution v with bounded 
initial data. The key point is that we get estimate fl84|) with precisely c°^{rj) 
in the exponent, rather than c°°{ri) + e, an issue that has become nontrivial 
because we lost the notion of self-adjointness when abandoning the Hilbert 
space setting. There is nothing to do in the case rj = rj^r, where we did not 
claim such an improvement. Also, the L'^{Ai) estimate for 77 = rjcr is clear 
because the operator is self-adjoint (with the Riemannian volume on the cigar 
as measure) in this case. 

We will distinguish the cases r] < rjcr and r] > rjcr- We rely on the fact that 
c°°{ri) > c°°{r]cr) and a trade-off between time decay and spatial growth that 
will be exhibited during the proof details. The key idea is that for s — )■ 00, 
the maximum principle 'almost' gives the correct growth rate c°°{ri) in the L°° 
norm. There is a remnant of the 0**^ order term that decays like 0(1/ cosh^ s) 
and that could give rise to secular terms; controlling this effect gives rise to the 
said trade-off between spatial growth and time decay. To control the estimates 
for small s, we use conjugacy and estimates with a better decay constant than 
c°°{r]) for data in the correct spectral space. 

The invariant projection corresponding to those eigenvalues that are 
strictly larger than c°°{ri) immediately gives rise to the decomposition f lS^ . In 
case c°°{ri) is itself among the L^m eigenvalues X^k, components in its eigenspace 
immediately satisfy (|8^ . So we can project this eigenvalue away, too, and need 
to prove only for v := Vres in the complementing space. Moreover, for 
the price of a small time shift (regularization estimate in Thm. 17. ip . we may 
assume Vq G C". 

Case 1: rj < rjcr- As remarked after f l79|) . multiplication by (coshs)^'^'' gives 
an isometry between L^(A^) and L^m(R"). Let us first conjugate the problem 
with (cosh s)'''"'"'' to obtain the self-adjoint operator L^^^ with L^(A^) initial 
data Vq := (cosh s)''~^'"''i;o- We find some c°°{ri) — e strictly above the next 
smaller of the L^(A^) eigenvalues (or, if there are no such eigenvalues, strictly 
above c°°{r]cr))- 

By the spectral properties of L^^^ in C"(A^), we conclude, uniformly in t, 

^ c||fo||c", this gives in particular the weighted estimate 
\v{t,x)\ < c(coshs)''--''e(^°°('')-")*||t;o||c- • 
So we choose 5 = e/ {rjcr — rj) > such that 
(85) suv\v{t,x)\ < ce^°°('')*||woi|c- ■ 

s<St 
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Flow of 

on -y G C" , 



or 



Flow of 
L„ on V 

'Icr 



L on (cosh s)''t; G = (cosh G C^_^^^ 

CL'^iM) ... since Tj <T^cr ■■■ C C"(A^) 

We get a spectral gap from the theory, 
but the estimates do not enforce sufficient decay at oo. 

A 

estimate applies to complement of eigenspace 

^^^ioT these eigenvalues 

^°°(^) K 

— possible further 1? eigenvalues 



^(% 



(in residual spectrum for C°~ 




In the gray wedge the spectral gap can compensate for the weight factor 

that is too weak. Outside, the maximum principle has to give growth 
control. 

Figure 3. Schematic outline of proof idea for rj < rjcr'- In the 
shaded cone, the spectral gap e compensates for the 'too slow' 
decay. In the exterior region, the estimate carries over by the 
maximum principle. 



We want to apply the maximum principle to obtain the same kind of estimate 
in the set s > 5t as well. To this end, v := e~^°°^'^^*v satisfies [dt — (L^ — 
= 0. While — c°^{rj) barely fails to satisfy a straightforward 
maximum principle - its 0**^ order term is positive according to (!25|) -, this 
term is < c/ cosh^ s < ce"^^*. So we let a{t) := c°°{T])t + ce"^*^* dt and define 
instead the new v := exp[—a(t)]v , which solves [dt — (L,, — a(t)))v = and 
satisfies a maximum principle. We get therefore, in the domain s > 6t: 

\v{t,x)\ = e'*W|y(t,a;)| < e"(*)c||t;o||L- < Ce^'^^^lt^oIlL- • 

Together we have shown ||'y(t)||L°° < C'e'^°°*^''-'*||'i;o||ca. With the regularization 
estimate \\v(t + l)||c" < C'||'u(t)||Loo from Thm. 17. the desired conclusion 
follows. 



Case 2: rj > rjcr- Here we let w^^* be the solution of an exterior initial 
boundary value problem that can be estimated with the maximum principle; 
then the remainder v — w^^* has compactly supported initial data, which allow 
for using the conjugacy again. This decomposition idea requires the use of a 
cutoff function x, which in turn entails writing some explicit projections to get 
back into the spectral spaces destroyed by the cutoff. Moreover, the remainder 
V — xv^^^ satisfies the differential equation with an inhomogeneity created by 
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the cutoff; however the effect of the inhomogeneity can be controlled, see ( 189|) - 
dnH) below. 

In detail, we first choose e such that c°°{r]) — 2e is still larger than c°°{ricr), and 
any larger pertinent eigenvalues buried in the essential spectrum for C° as 
well, and we choose 6 such that 6{ri — rjcr) = £■ We then we choose (ti,Xi) 
(with the geodesic radial coordinate si) and let v'^^^ solve {dt — L^)-?;^^* = in 
s > Si + S{ti — t), < t < ti, with boundary data v^^^ = if s = Si + —t), 
and initial data Vq^^ = vq for s > si + ^ti + 1, and smoothly cut off to for 
s G [si + 5ti,si + 5ti + 1]. 



-ext 



{dt - L^)t;^^* = 



V ^ cutoff Vq s 

< r> 

Si + dti 

To prove the existence and uniqueness of the solution to this exterior boundary 
value problem, we first map the exterior of the cone to the exterior of a cylinder 
by means of the map s + (5t — (si + 6ti) = a, with a the new radial coordinate, 
and then refer to Cor. 16.21 The estimates obtained there are uniform in (ti, si) 
as long as Si is bounded away from 0. Indeed, with v{s,ijj,t) = V{a,u,t), the 
equation {dt — L^)f; = transforms into {dt — 6da — Lr^)^ = 0, and is the 
same as in f l25|) . except with replacing dg. 

With a{t) := c°°(r/)t + c /J e-2(«i+5(*i-*))rft, the function := e-'^^t;'^^* sat- 
isfies a maximum principle again, and we get 

(86) |r^*(t,a;)| = e"(*)|F^*(t,a;)| < e^'^'^lvoh^ < Ce^'^^^it^ollL- 

with a constant C = C{6) that is easily calculated to be uniform in (ti,Si) 
as long as si is bounded away from 0. In particular, ||'y'^^*(t — < 
(7gC°°('7)(<-i) Having been constructed by means of a cutoff function, 

^ext j^g^y j^Q^ ^YiQ appropriate spectral space. But with a regularization 

estimate and successive projection, we can conclude that 

(87) ||P,r^*(ti)||c.. < C\\v'-\h)\\c^ < Ce^'^^'^^'^WvoWc'^ . 

Now take a cutoff function x ^ C°°(R) for which x(cr) = if a < 1 and 
X{a) = 1 if (T > 2. Let 

w»* ■=v- P^(r^*x(s - si - 6{h - t))) = P^{v - xv^""') , 

where we have used in the last step that v is assumed to be in the range of 
P^. The initial data of ■u*"* := v — x-u^^* are compactly supported. So we can 
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conjugate, letting 

(88) := (coslis)^-^" t;'"* and w"^' = P^.^w;'"* = (coshs)'^-^-t;'"* . 

We have 

, . =-(coshs)^-^-(9i-L,)(xt;^-*) 
^ ' = -(cosh s)'^-''^^^ [dt - L^, x]^^'^* 

= (cosh s)''-''- (x(a)t;^^* + 2x'i(T)dsV^^^) =: f + d^g =: h 

where 

g := 2 (cosh and 

/ := (coshs)^-^='-x(V)t;"^* - 2ds{x'{(r) (cosh s)"""") v'""' 

and, using m and m, X = - fe°^(r/) tanh + x" e C'o~]l,2[, and 

(T := s — Si — 5(ti — t). 

We'd like to argue that 

~v^-\t) = S,„,(t)^-*(0) + £ S,„,(t - dr , 

5-*(t) = S,„,(t)S-*(0) + S,„,(t - r)P,^>(r) dr . 

However, we have proved insufficient regularity for /i(r) to be an initial value 
for the semigroup estimate. While this could be fixed by stronger regularity 
estimates, we can more easily average over estimates for initial times to ^ [0, 1] 
and use Thm 17.1] directly. At a formal level the argument is written as follows: 
for t > 1 and all to ^ [0, 1], "we have 

(90) ^-*(t) = Sr,Jt - toW'^'ito) + f S,„,(i -to- r)h{to + r) dr ; 

Jo 

so by averaging over t^, we have 

i,'^\t) = J^S,Jt~to)i,'^\to)dto 

+ lo lo'^ ^vcr it -to- 'r)h{to + t) dr dto 
+ lo It-i" ^v.r it -to- T)h{to + r) dr dto 
= S,^,Xt - 1) Sr,M - toW'^'ito) dto 

(91) + Sr,At - 1 - r) S,,,(l - to)M^ + ^o) dto) dr 
+ (/o S,,,(l - to - r'Mt - (1 - to) + r') dr') dto 

=: S,,,(t-l)yo 

+ C Sr,Jt - 1 - r)V,{r) dr 
+ JoV2{l-to)dto. 

The fact that the integrals in flHT]) still manifestly apply the semigroup to a 
potentially singular distribution h is of no concern for the argument. The Vo,i,2 
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terms are simply solutions of the inhomogeneous equation with homogeneous 
initial data, and they can be estimated by Theorem 17.11 

Vo is the solution at time 1 of {dt — L^^^)V^ = -y*"* with initial data 0, and so 
we have, using ( 15^ . (15^ along with Thm. I7.H 

||K)||c"(>!) < C'||'^'™*IU°°([0,l]xA4) ; 

||'«^*"*||c"([0,l]xX) < C'(||'l'*"*(0)||c'^(X) + II/, fi'||L°=([0.1]x>!)) 

< Ce(^i+'5*i)('?-^-)(||^o||c"(A,) + \\v'^'\\L^ao,l]xM)) 

< Ce(^i+^*i)(''-''-)||i;o||c-(>() • 

Vi(r) is the solution at time 1 of {dt — L^^^)V^ = h{T + ■) with initial data 0. 
We have 

ll^l(''")llc«(A^) < C||/, 5'||L°=([r,r+l]xA^) 

^2(1 — to) is the solution at time 1 — to < 1 of ((9f — L^„.)^ = h{t — (1 — to) + ") 
with initial data 0. We have 

||V'2(1 - to)llc-(AI) < C\\f,g\\L'^(^[t-(l-to),t]xM) 

With these estimates established, we can now apply the spectral projection 
and conclude 

Jo Jo 

hence 

||^^"*(t)||c"(A^)<C7e*(^°°('')-'^)||K,||c"(>i) 

+ C/o"'e(*-i-)(^-(^)-2^)||\4(r)||c.(A,) 

+ <^SUPigg[0,i] 11^2(1 -to)||c-(A4) • 

The first and third term can immediately be estimated as C exp[t(c°°(?7) — 2e) + 

{si+6ti){r]-r]cr)]\\vo\\c'-{M) iind C exp[{si+6(ti-t)){T]-r]cr)+c'^{r])t]\ML°°{M) 
respectively, and for t = ti, they are controlled (with one e to spare in the first 
term) by C exp[tic°° (rj) + si(?7 — r/cr)] ||^o||c°(ai)- The second term is estimated 
similarly to be 

<C [ exp[{c°^{r])-2e){t-l-T) + {si+S{h-T)){r]~r],r)+c°°iv)r] rfr||t;o||c^(A,) , 
Jo 

which, for t = ti, is < C exp[tic°°{ri)+si{ri—ricr)] Jq^ ^ ^-siti-i-r) ||f;o||c''»(_A^), 
and the integral of order 1/e can be absorbed in the constant. 

So we have proved 

||^^"*(ti)||c<^(A4) <Ce*^^°°('')+^i(''-''-)||t;o||c"(A,) ; 
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and in particular \ff''\ti,Xi)\ < Ce^'^^^'^^+'^'^'T'^^^l With (EH]), this implies 
\v^"^{ti, xi)\ < Ce^^'^°°^'^\ The constant is uniform in as long as si 

is bounded away from 0, say Si > 1. But for s < 1, the estimate with 
Si = 1 still gives a;)| < Ce^'^^°°^'^\ We review (186 p . which guaranteed 

ll'^'^^*(^i)IU°° — II-UoIIlo", and since the projection operator is bounded 

from into itself, we also have ||(P^(xt)'''*))(ti)||L- < Ce^'^^'^^MI^olU-- In 
view of v{ti, Xi) = Xi) + (P^(x'y'^'^*))(ti, a^i), which we have from (190|1 . 

we conclude 

\vih,x,)\<Ce'^'^^^^'^\\vo\\c^ 

uniformly in {xi,ti). This establishes ||'y(t)||Loo < Ce^'^^^^^\\vo\\c'^; and with 
another invocation of regularization, we can get the same estimate for the C" 
norm on the left hand side. □ 



9. Proof of Theorem 11.11 

Let po be an initial density which satisfies the assumptions of Theorem 11.11 
After scaling we may assume that its mass agrees with that of Ub{x) for B = 1. 
We shall actually prove the asserted decay (fTOj) in a norm || ■ ||c"(X) stronger 
than II • ||L°o(j^n). As long as m > mi = the center of mass of po is defined. 
We may shift po, such as to bring its center of mass to 0. We then claim that 
the statement of the theorem holds in particular by choosing z = 0. The 
hypothesis m > mi is automatically satisfied for n = 1, and also whenever 
m > rrin- In the other cases, we either cannot define a center of mass or else 
will not bother to shift it, and we claim that in those cases, the theorem holds 
with any z. For quick reference, here are the relevant eigenvalues of L in C" 
(unweighted) in the case n > 2: 

Xqq mass 




Aoo = 

Xiotransl eigenvalues: Aio = -(p + n) 

Xqi scaling Aqi = —2p 

In the case n = 1 (where m„ = 0), the picture looks like this: 
?inn mass 

— . , Aoo = 

eigenvalues \ t , i\ 

^transl Jen n = 1: ^- = -(^ + 1) 

Aqi scalmg Aqi — — 

In either case, the onset of the {£ = layer of the) essential spectrum is right 
at Aoi- 



The evolution conserves mass and center of mass. In the variables of eqn. fl2T|) . 
these conservation laws amount to J{v — 1)ub dx = and j{v — 1)ub xdx = 0. 
In other words, using the inner product of Lf^m makes v — 1 orthogonal to 
the eigenfunction t>oo = u]f"^ and, if m > mi, to i^io = ru]f"^Yi^{u) for 
fi = 1, . . . ,n. 
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By ([3]), the statement of the theorem, with z = 0, is equivalent to 



sup e 



-Aoitl 



v{x, t) — 1| < OO 



recalhng Aqi = —2p. Let us rewrite fl2T|) - fl22|) as 

(92) 



idt-L)iv 



f 2(n-l) d , , , 2. 



m 



+ n 



+ 



2(m + 1^ 



tanh s— V 

os\m 



m] 



m 



n- 



1 + (cosh s) {mv — f + 1 — m) . 



With the trivial substitution w 
form 



1, its structure becomes clearer in the 



(93) {dt - L)w = Uf{w)w) 

where f{w) = ((1 + w)"^ — 1 — mw)/'mw is analytic and vanishes at the 
origin, and where L = L — Y^(tanh sdg + n — tanh^ s) can be read off 
the right side of f l92|) . We can apply Thm. 17.81 Let us write Wj := w{jT) G 
C"(A^), for a constant T chosen sufficiently small for the contraction estimates 
of Chapter [TJ specifically from Lemma 1731 to apply. Also let S(t) := exp[tL] 
be the semigroup generated by the linear operator L. We write 

(94) Wj+i = S{T)wj + g{wj) , 
where 

C"(M) g{w) e C'iM) 

is smooth and vanishes at the origin together with its derivative. More specif- 
ically, since the map Wj Wj+i =: F{wj) is smooth according to Thm. 17. 8[ 
and F(0) = 0, DF{0) = S(T), we can write wj+i = £F{awj) da = 

S{T)wj + [jQ{DF{aWj) — S{T))da]wj =: S{T)wj + G{wj)wj. Here G is smooth 
with values G{w) being bounded linear maps C", and G{0) = 0. This 

gives the claim. 

In the case m > rrin, let P be the S-invariant projection to the L^m(R"') 
complement of the eigenspaces for Aqo = and \iq = —{p + n). In the case 
m < rrin, let P be the S-invariant projection in C"(A^) to the complement 
of the eigenspace for Aqo = 0. It is well-defined according to Theorem 18.51 
We have seen that the Wj are in the range of P, and therefore ||S(t)wj|| < 
(^gAoit||^^.|| Tffm^ a constant C independent of t, where || ■ || refers to || ■ ||c"{A^)- 

We introduce a small quantity r < 1 and assume ||wo|| is less than r divided by 
a constant supplied from Lemma 17^ This lemma then ensures that ||wj||L°° < 
r for all r. Lemma [7.51 applies to our operators L, L and nonlinearity f{w)w, 
and it controls [[(^(w)!!^^ < -ft'llwUcQ llw^lU"" for IIwUloo < r. Now K stays 

62 



fixed, while we still may impose further smallness requirements on r. We 
denote by A < whatever spectral gap applies to the parameter m under 
consideration. Iterating gives 



(95) 



Here, C is the constant obtained from Thm. 18. 5[ We know that 

II I! ^ • r I! I! 1 ^ 1/21! 111/2 

||Wj_i||Loo < mmjr, < r ' 

and will prove inductively that < 2Ce'^'^"^||wo||c"- The start being 

trivial, the induction step in fl95|) follows if we can ascertain that 

1 + ^;^g-AiT^l/2(2C<)3/2e3A0-l)T/2 < 3 . 

i=i 

By convergence of the geometric series e^^'^^'^, this requirement can indeed 
be met by choosing r sufficiently small. 

Apart from the adjustment of P, the proof works the same in the full range 
of m (with A + e in case m = 1712, per Thm. 18. 5p . 

Finally, let's deal with the case m = m2 under the moment hypothesis (|9]), 
which, in view of ( 128|) and (130|) . amounts to fj^wo{yY dfi < 00. (Recall dii 
denotes the Riemannian volume element on A^.) We will argue shortly that 

(96) \\w{t)\\L2(M) < Cc^^WwoWl^m) ■ 

We also note that ||S(T)t(;||cc(_A4) < C||w||l2(_a4) by regularization: indeed, 
this is a regularization estimate for a uniformly parabolic linear PDE in each 
coordinate chart Ui separately, with the ||w||i2([;j) all being trivially estimated 
by the global ||w||j;^2(^); recall that the distortion between the metric on the 
Ui and the euclidean metric is uniformly bounded. 

From these two ingredients we conclude, in view of the close-to-sharp 0(e'^^"''^^*) 
decay already established, that 

\\wit + T)\\c^ < ||S(T)«;(t)||c;<. + l|^7(«;(t))||c« <C7e^*lkolU2(X)+Ce2(^+^)*, 

and the claim follows since 2(A + £:) < A. So all we have to establish is fl96]) . 
and we do it by direct calculation from the PDE. 

Letting h{w) = ^^^""J""^ = w + 0{w'^) and H{w) = ^w^ + ©(w^) its anti- 
derivative, we use E{t) := Jj^H{w{t)) dfi as a proxy for the L^-norm. For 
the technical reason of properly caring about a flux term at infinity, we will 
also need E{r,t) := f^^^^ H{w(t)) dfj,, with integration over the euclidean r- 

ball around the origin (equivalently the geodesic arsinh r-ball) . Formally, we 
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calculate 

dtE{t) = I h{w)wtdii 



M 



(97) = I h{w)Lh{w)dn+ I h{w){x-V){w -h{w))dn + 

J M I — m J 

\x\^'^ {w — h{w)) h{w) djj, . 



M — III J 

2 f ( 2m^-™ 

+1 / r^-T^ ' 

^-tn Jm^ 1-m 

Returning to the euclidean volume element d^x = (1 + |a;p)"/^(i/x, we want to 
integrate the second term of this sum by parts to avoid derivatives of w. With 
h{w)'V{w — h{w)) — 'V{H{w) — we obtain for n 7^ 2 (with routine 

modifications in case n — 2) 

dtE{r,t) = I h{w)wt djji = I h{w)'Lh{w) d/j, + 

JB(r) JB{r) 

2 



+ 



(l-m)(2-n) Jb 
2 



/ V ■ Uh{w) - \h{wf) V(l + d^'x 

JB(r) ^ ' 

f {H{w) - lh{wf) A(l + |a;n-"/2+i ^x 
Jem 



1 -m)(2-n) 

/ in — \x\'^ ) (w - h(w)) h(w) d/j, . 

JB(r) ^ l-m / 



1 — m 



The last integrand is estimated as 0{w'^), or 0{wH{w)) = e'^^'^^^^0{H{w)). 
The second last term is integrable even over Ai, or can be estimated in the 
same manner as the last term. We obtain (with dS the Riemannian surface 
element) 



dtE{r, t)< ! hLh dix + Cie(^+^)* / \w{t)\^ dS + C2e^^"''^*E{r, t) . 

JB{r) JdB{r) 

Letting -f{t) := - /J C2e^^+=)* dt, we get 

dt{e^^*^E{r, t)) < e^^*) / h{w{t))Lh{w{t)) d^ + 



(98) •'^('^^ . 

+Cie^(W+^)* / \w{t)\^dS 

JdB(r) 



and therefore 



e^(*)£;(r, t) - E{r, 0) < [ [ h{w{t))Lh{w{t)) dn dt + C\\w{t) \\loo . 

Jo JB{r) 



Now from ||wo||i2(_yv/() < cxd, we get £'(r, 0) — )■ £"(0) < cxd as r — )> cxd. If 
w{t) G L^, we conclude ^^^^^ h{w{t))Lh{w(t)) djj, — )> h{w(t))'Lh{w(t)) d/i < 
Xh{w{t))'^ dfi < 0. Without assuming w{t) G L^, we still maintain the in- 
equality limsupj.^^ J^^^^ h{w{t))Lh{w{t)) d/i < by approximation, and this 
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guarantees the finiteness of liiiir^oo e'*'*-*''-E(r, t), and therefore w{t) G L^(A^) 
for every t. 

Knowing this, we can actually justify the formal calculation (P7|) and proceed 
with the same estimates and integration by parts done for E{r,t) to obtain 

dtEit) < \ [ h{wit))Lh{wit)) d^i + Ce^^+'^'E{t) < (2A + Ce^^+'^')E{t) . 
Jm 

The estimate E{t) < Ce^^*£'(0) follows from this; and then immediately the 
same estimate for ||w||^2(_yvi)) hence ( p6|) . 

10. Asymptotic estimates in weighted spaces: The case m < 

In the previous section, we have seen in particular that for m < 1712 = we 
get a convergence rate of 0(l/r), and no better, to the appropriately shifted 
Barenblatt with respect to the relative L°° norm. This is the same rate as was 
obtained by Kim and McCann |3T] via Newton potentials, and later extended 
to a wider range of m by the quintet [7] at the cost of an extra e > in the rate 
exponent. The slightly larger spectral gap in the weighted Hilbert space of [20] 
(where Aqi was not an eigenvalue any more for m < m2 and the continuous 
spectrum started further below) was not actually expected to be dynamically 
effective from that setting, because the linearization formalism relied on the 
existence of second moments, in contrast to the reality in the case m < 1712. 

In the present setting, the absence of a weight in has changed the spectrum 
in comparison to |20]. To recover that spectrum, we need to study L in C^^^, 
or equivalently, the operator L^^^ in C". Since rjcr = ^ — ^ < for m < 1712, this 
means we need to study the case where initial data deviate from Barenblatt less 
than what would automatically be achieved by Vazquez' result dS]). For initial 
data with this special tail behavior, we can indeed get improved convergence 
rates from a study of the conjugated operator. These better convergence rates 
will apply to appropriately weighted norms. 

We are careful to point out that the nonlinear evolution was constructed for 
unweighted norms only: this is in accordance with the fact that, in the relative 
L°° setting, unweighted norms keep a uniform distance from the singularity 
u = of the Fast Diffusion Equation ([1]) and allow, on a technical level, 
uniformly parabolic estimates on the whole space. (However, in the present 
case m < m2, Thm. 17.81 does carry over, as noted near the end of Section [3) 
In contrast, the linearized flow can be studied with the same ease in weighted 
and unweighted spaces. 

Let's have a quick look at the spectrum of L^^^ in this range of parameters, 
where the only eigenvalues are Aqo = and Aio = —{p + n), and the onset of 
the essential spectrum is Xq°^^ = — (| + 1)^: 
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In the case n = 1, one has mi = 0, so Barenblatt always has 1st moments and 
we can shift the center of mass. We also need to do this because Aio is above 
the onset of the essential spectrum. 




eigenvalues of "^^^ _^ , . 

L when n - 1- ^° " + 

^Vcr wnen il — i. ^^^^^ _ ^ 

In the case n > 2 the center of mass may or may not be defined (depending 
on whether m > mi or m < mi), but there is no need to shift the center of 
mass, because the intersection of — (| + 1)^ with Aio = —{p + n) occurs at p = 
2y/n — 1 > 2, hence Aio is below the onset of the essential spectrum (ranging 
from —1 to —4 monotonically) in the whole parameter range mo < m < m2. 

As we have resolved to work in the space C^^^, we let w = {cosh.s)~^''''w, and 
now equation fl93|) becomes 

(99) (dt - L^Jw = L,^, (/(I + (cosh sy^^w)w) 

Since rjcr < 0, we can estimate || (cosh s)'''=''w||ca(A^) < c||w||c"(a^), and this is 
why Theorem 17.81 carries over to this case along with its proof, as remarked 
at f l77|) . We now can repeat the proof from Section |9] almost verbatim. We 
conclude: If the initial data wq satisfy 

(100) (cosh s)-''^'- Wo e L'^iM) , 
then 

(101) supe(-^o°-"--)*||(coshs)-^^'-u;(t)|Uoo < oo 
t>i 

for every e > 0. (The same statement works also for the C"-norm.) Here 
-A§°^* = (1 + 1)^ = (t^ + 1 - i)^- If we put nL^ in the hypothesis ffTOOj) . 
we get conclusion ( llOip without the e. 

We therefore get the following: 

Theorem 10.1 (Fine Asymptotics for m < m2, with restricted initial data). 

Assume n > 2 and m G ]mo,m2[ = ]^^);^[; cr else n = 1 and < m < 
m2 = |. Further assume that the mass of po is one, and if n = 1 also that 
the center of mass of po is 0. Let S := = y — (note 5 < 0). If 

Ps(0, ■)'^|1 - po/Ps(0, ■)! ^ L°°, then, for every e > 



(102) limsupsupr 



7-e 



y 



pB[r,y) 
PB{r,0) 



pB[r, y) 



< oo 



where 7 = ^(| + 1)^. Likewise, we get the same conclusion without the e, 

provided we add the hypothesis that ps(0,-)''|l — po/Ps(0, ■)! ^ ^^so 
square integrable with respect to the measure (1 + \x\'^)~''^^'^ d"'x. 
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A graph for 7, 6 in dependence on m can be found in Figure [H The tremen- 
dous improvement over the 0(l/r) convergence rate estabhshed by Kim and 
McCann [3T] and in our present Thm. 11.11 depends on the strong decay as- 
sumptions for the relative deviation of the initial data. 

Also note that our Theorem 110.11 is a generalization of the comparison esti- 
mates by Vazquez (of which we use a modified form as Lemma [7.6p . Namely 
this comparison asserted that if a solution is close to Barenblatt in relative 
L°° norm, then it stays close to Barenblatt in this norm. Here we have shown 
the same in a stronger norm, showing that proximity to Barenblatt in a norm 
that enforces an even closer fit with Barenblatt in the tails, is also preserved 
under the dynamics. 



11. Higher asymptotics in weighted spaces: The case m > 

n+2 

Proof of Theorem 11.21 and its corollaries. 

This section is devoted to proving Theorems 11.21 and lll.li When applicable, 
the latter extends the conclusions of the former to decay rates A outside the 
range ]2Aoi,Aoi]. For an example of its apphcability, see Remark 11.91 The 
restriction A > A + Aqi appearing in this refinement arises from a limitation 
of our method, which relies on a weighted control of the quadratic term in the 
nonlinearity by a product of a weighted norm (I103P decaying at rate A and an 
unweighted norm decaying at the rate Aqi (given by Theorem ll.ip . 

It seems convenient to repeat Theorem 11.21 here: 



Theorem 1.2 (Higher-order asymptotics in weighted Holder spaces). Fix 

p = 2(1 — m)^^ — n > 2 (equivalently m G ]m2, \[) and A € [Aq°"*, Aqi] = 
[— (| + 1)'^,— 2p] subject to the condition 2Aoi < A. Ifu{t,x) is a solution 
to 11]) with center of mass and linit_^oo \\u{t,x)/uB{x) — = 

both vanishing, then there exist a sequence of polynomials {uik{x))ik, each 
element of which either vanishes or has degree £ + 2k + such that 

[B + \x\^){u{t, x)/us{x) - 1) - Ea<a,,<o Uf,k{x)e^' 
(B+|a3|2)(''+'-V(P+')'+^4/^ 

as t ^ 00, where the sum is over non-negative integers fc,£ G N for which 
Xik defined by (llip lies in the interval ]A, 0[, and for which £ < 1 if n = 1. 
The functions uik{x)/{B + \x\'^) lie in the Xit eigenspace of the linear 
operator on C^^^(7M), and the norm \\ ■ WciM) — II ' IIl°°(R") '-s 

defined by p3p . 

Theorem 11.1 (Higher-order asymptotics in weighted Holder spaces). Fix 

p = 2(1 — m)^^ — n > 2 and A G [Aq°°*, Aqi], and choose A < Aqi such that 
A + Aoi < A. If u(t,x) satisfies the hypotheses of Theorem and either 
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0{e^') 



A = Aoi or else we have the stronger hypothesis 



(103) limsupe 



-\t 



t—>ca 



u(t,x)/uB{x) 



< oo , 



then the conclusion f lT4|) of Theorem remains true and u^k = for Xek > A. 

Proof of Thm. In case A = Aqi = — 2p, the weights cancel in f lT^ . and 
the result was established in Section |9] during the proof of Theorem 11.11 The 
improvement to A < Aqi is based on the same principle as Section [101 except 
this time rjcr > 0, and the situation is dual to the previous one: we have 
a chance to obtain faster decay by introducing a weight which relaxes the 
strength of the norm as s — )■ oo. We need not conjugate all the way to rjcr, 
but instead have the flexibility to choose t] G [0, to balance the error term 
in Theorem 11.21 against the severity of this relaxation. Indeed, the choice 

(104) V = Vcr- VA - A-°* 

makes the radius Aq°^* of the essential spectrum of given by Theorem (18.20 
coincide with A G [Aq™*, Aqi]. Choosing rj > rjcr would only prove the statement 
with a weaker norm, without giving better rates. 

Conjugating with rj > however requires a reconsideration of the proof of 
Theorem I7.8[ because a priori, uniform parabolicity cannot be expected for 
weights that allow growth relative to the Barenblatt. The redeeming feature 
will first be that we do retain the unweighted relative hypothesis on the 
initial data; then since the nonlinearity enters through a term f{w)w, after 
conjugacy w = {cosh s)~'^w, this term is still f{w)w, and the unweighted 
estimates continue to control the nonlinearity, whereas linear theory applies 
to the weighted estimates. We now carry this out in detail: 

We again study equation ( 193|) for initial data wq G (in particular bounded), 
but using weighted C° norms with the more permissive weight r] > from (11041) . 
Letting w := (coshs)~^w, Eq. (1^ becomes 

(105) {dt-L^)w=L,{fiw)w) 

with = — (tanh s dg + n + {i] — ^^^^"^ ^) with the same f{w) 
as before. We can use the fact, from the unweighted norm, that a solution 
w G C"(A1t) to ( l93l) exists; in particular the conjugated w{-,t) will also be in 
CiMr) still (even in C°^). 

We obtain an analog to Lemma 17.51 but stated specifically for the FDE on 
the cigar manifold, because the a priori estimate on relative norms from 
Lemma [7.61 is used: namely, we claim 

Lemma 11.2 (Linear approximation of nonlinear semiflow; weighted norm). 

Let w solve the homogeneous linear equation ( 15^ . namely (9t — L^)w = 
for initial data Wq, where is given by (1251) and can be written, in local 



coordinates, as L^zi) = dfji^a^^w) + diiUw) + cw. Let w solve the quasilinear 
equation fllOSp . for the same initial data Wq, and w = {cosh s)'^w. In local 
coordinates, we can write L = L + 9^ o 6* + c, and f is a smooth function from 
an interval about into R satisfying /(O) = 0. The coefficients are smooth. 

Then, for sufficiently short time T, there exists a constant K (uniform as 
T — )■ j such that we have the estimate 

(106) \\w - w||c«(A^t) ^ K\\M\l^{Mt)\\M\c'-(Mt) 
and from it the time-step estimate (with a different K): 

(107) ll^'(^) - w{T)\\c'^(^M) < KWwoWc-iM) ||w^o||l°°(>!) ■ 

Proof of the Lemma. The proof is modeled right after the proof of Lemma [7751 
We can basically copy equations flTTl) and (1721) almost verbatim, the only 
changes being that all w^s have tildes now, except the ones inside /(■); and the 
insignificant fact that coefficients implicitly depend on the conjugacy param- 
eter f] (not to be confused with the partition of unity t]i). The proof ensues 

as before, with the estimate \\fiw)w\\c'^(MT) ^ -^II/('^)IU°°(a^t) II'^IIc'"(A4t) 
provided by Corollary I5.9[ □ 

Proof of Thm. li.i^ continued: As before in Ch. [HI we write Wj := w{jT) G 
C"(A^), where T is sufficiently small for Lemma [11.21 to apply; and we write 
S^(t) := exp(tL^) for the semigroup. We continue to use Thm. 17.81 for the 
unconjugated flow and copy from flM|) and the paragraph following it that 

Wj^i = S(T)wj + G{wj)wj . 

After conjugating, this becomes 

(108) Wj+i = S^{T)wj + Gr,{wj)wj , 

where Grj is a smooth function on C" with values Gr^{w) = (cosh s)~^ o G{w) o 
(coshs)'' being bounded linear maps on C".^(A^). The estimate fl70l) guaran- 
tees that G>n{wj) is also a bounded linear map on C"(A^) (namely with norm 

< fr||tyj||2,oo). 

From Theorem 18.21 the only spectrum of closer to zero than Ag°^* = A 
consist of eigenvalues ffTTl) of flnite multiplicity and indexed by non-negative 
integers £, /c G N such that £ + 2k <p/2-\-l — Irj — r/crl- Enumerate the 
eigenvalues X^k which lie in the range ]A, 0] by A^ < A;_i < • ■ ■ < Ai, counting 
them with multiplicity. Set A^+i = A by convention, whether or not this is an 
eigenvalue. 

Before continuing the analysis of fill 01) . a brief remark about the heuristics of 
the proof strategy seems useful: We are not constructing a 'slow' or 'pseudo- 
center' manifold corresponding to the eigenvalues Aj using invariant manifold 
methods with respect to which one would naturally measure convergence rates 
of the transversal 'fast' dynamics. Rather we assess convergence rates in terms 
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of distances to the linear eigenspaces using spectral projections. This is appro- 
priate since the spectral ratio A/Aqi is less than 2, the order of the nonlinearity. 
For larger spectral ratios, the deviation of the slow manifold from its tangential 
space, rather than the eigenvalues alone, would indeed dominate the asymp- 
totic behavior of the 'fast' dynamics transversal to the slow manifold. This 
is also the reason why products of eigenfunctions and linear combinations of 
eigenvalues do not occur in the statement of the theorem. (Fig. H] explains this 
heuristics in a simple 2D model.) 



\\ t 


X = \x 

nonlin. system: . . , 9 

y = Ay + 

/ A. ^ X = Ax: x(t) = xne'^* 
.. ■ / hnear system: . . )/. 

A slow mfd for the linear system 




eigenspace for A, A slow mfd for the nonlin. system 



Figure 4. A simple 2D model illustrating the heuristics of 
working with linear eigenspaces with A < X < in Thm. \1.2\ 
rather than slow manifolds. Linear asymptotics is XqC^^ + 0(e^*) 
(solid curve). Nonlinear asymptotics (dotted curve) should in- 
volve e^^* etc., which gets absorbed into 0(e^*) if A/A < 2. 



Returning to (11081) . let Q*^*-* be the spectral projection onto the one-dimensional 
eigenspace corresponding to the eigenvalue A, of L. Take Q = J2i=i Q^^\ 
let P be the complementary projection. 

Setting Q^,*^ = (cosh s)~^ o Q*^'^ o (cosh s)'' and = (cosh s)~'^ o P o (cosh s)'^, 
we obtain 

P>j + 1 = Sr,{T){Pr,Wj) + (P^^ O G{Wj))Wj . 

By iteration, we obtain inductively 

fc-i 

j=0 

We set -di+i := exp[A;+iT], so ||S^(T)'=|| < C^f^^ on the range of P^. We 
also use what we have already established in Thm. II. 1^ namely that Hu^j < 
C-d^ with -d := exp(AoiT), with a similar estimate following for wj. Then we 
conclude 

k-l 

(109) IIP.^fcllc'^cA,) < c^t,{\\p,wo\\cnM) + ■ 

j=0 
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Since ^'^/^i+i < 1 (here we use 2Aoi < A), the parenthesis is bounded inde- 
pendently of k. So we get 

(110) ||P>fc||c-(A^) < C2^9f+i . 

A similar estimate can be made for the complementary projection 

i=l 

On the Aj eigenspace, S^(T) operates as the scalar {}i = exp[AjT], hence 

k-l 
j=0 

leads to a 2-sided Cauchy sequence estimate 



k+q 



9-1 



We conclude Vi := limfc_>.oo Qlj Wk/'&i exists, and Vi must be an eigenvector of 
in the Aj eigenspace. The convergence rate follows from this estimate: 



(111) 



,~. „Qfc I 



< C3{^ 



2\k 



- (coshs)''fj, we use from 
Ui{x)/{B+\x\^) 



Returning to the unconjugated functions, with Vi 

Theorems 18. II and 18.21 that this eigenvector Vi{x) = Ui{x)/{1:5 + \x\'') is given 
by a degree i + 2k' polynomial Ui{x), and Aj = Xik'- Since Ub{x) has the same 
mass and the same center of mass as the initial data by hypothesis, we have 
uoq{x) = = uio{x) as remarked before equation fl92l) . hence ^ -|- 2fc > 1. 

The last two inequalities combine with the identity 1 = P + E!=i Q''*'' yield 





X 




- E Li 


{B + 


X 


2)l+r?/2 



(cosh s)"^ 



Our choices r/ = r/^r - V A - Afj""* = | - 1 - v^A + (| + 1)^ and ^9j = exp[Ajr] 
for T sufficiently large convert this into the desired estimate f|T4|) . This con- 
cludes the proof of Thm. 11.21 □ 

Proof of Thm. \11.1\ Since the case A = Aqi is Thm. 11.21 we assume A < Aqi 
and hypothesis f llOSp . The proof of Thm. 11.21 carries over with only minor 
changes: 

We keep conjugating with the r] from fll04p and instead of using \\wj\\ < 
C exp[AoiT_7] =: 'd^ obtained from Thm ll.l( we use the stronger f llOSp . \\wj\\ < 
Cexp[ATj] =: -dK As in the previous proof, we obtain f llOQp . only with 

{'d'dy instead of ■i?^-', and use < 1 from A + Aqi < A. The same 

change applies to the estimates of Q. 
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It is easy to see that the largest eigenvalue A^^ for which the polynomial 
Uik{x) 7^ is non- vanishing cannot exceed A without ffT^ contradicting ei- 
ther (11031) (in case A < Aqi) or Theorem 11.11 (in case A = Aqi). □ 

Corollary 11.3 (Coefficient formulae). The coefficients appearing in the as- 
ymptotic expansion (IT^ from Theorems and \ll.l\ are given by Uik{x) = 
i^ik{\x\) J^fj.'^^kfiyefiix/lxl) where ipek andYi^ are the polynomial eig en functions 
of H and spherical harmonics from Theorem \8.1\ and 

(112) Qfc^ = lim e-^"*(M(t, ■) - Ms,^£fc>'^M)L2(Rn)/||V^^fcl"£M«'""^'"' 



B llL2(R") • 



For k = ^ i, the ipeoYin are the homogeneous harmonic polynomials of degree 
i and the expression under the limit flll2p is independent of time, hence 



(113) ueo{x) = \x\'J2'^i^^{- 



X 



x\ 



{u{0, •), V '^oy£^)L2(R") 



l-m/2 ||2 



Proof. Rewrite f|T4|) in the form 

where we have L°° control on the factor 0(e^*), and Uik lies in the \ik eigenspace 
of the self-adjoint operator H on L^a-m spanned by the family of degree 

Ug 

i + 2k < p/2 + 1 polynomials {;ip(.kYe.n)^i- Multiplying this expression by 
g-Aftt^^^Y^^ and by = {B + |a;|2)-(n+P+2)/2^ integration yields 

e-^"*(M(t, ■) - Ms,^£fc>"^^)L2{R") = iu,k,'4^ikYt;)L\ + 0(e(^-^^'=)*) . 

The remainder term vanishes in the limit t — t- oo to establish f lll2p . 

In case = 7^ £ we have '?/'£o(|^k|) = \x\^, hence {'iptoY^^}^ are the homo- 
geneous harmonic polynomials of degree I. These integrate to zero against 
the radial distribution Ub- It remains only to show time independence of the 
expression under the limit flll2p to complete the proof of flll3p . Transforming 
back to the original variables ([5]), we see that the integral of p(r, y) against any 
harmonic polynomial of degree less than p is independent of r from the evo- 
lution equation ([1]); the spatial decay |Vp(r, y)| = 0{l/\y\'^'^'P'^^) in e.g. [311 



Corollary 9] justifies the integration by parts. Thus 
j p{T,y)\y\'Y,,{^^dy = j w(ln(l + 2pr)'/^ ^)l(l + 2pr)^a;|%(^)rfa; 

is independent of t = ^ln(l + 2pr). Since A^o = —2pf3i from ([3]) and (fTTj) . 
this establishes (11131) . □ 
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Proof of Corollary | i. 7[ Let us quickly review the top of the spectrum from 
Theorem l8.2l (omitting the eigenvalues rendered ineffective by the mass / center 
of mass adjustments): 

Aoi = -2((1 - m)-^ -n) = -2p 
effective eigenvalues A02 = — 8((1 — m)"^ — 1 — n/2) = —4p + 8 
of hr, when n > 3: A20 = — 2p — 2n 

^cont ^ + 2)2/4 

The same picture, with the interval [m6,m„+4] shrinking to a point, applies to 
n = 2, whereas in the case n = 1 we have the intersection of Aq""^* and A20 at 
p = 2 + 2^2. 

Apart from Aqi = —2p, which forms the top of the spectrum of L^^^ in the 
full range p > 2, the spectral gap is given by A = maxjAg""^*, A02, A20} = 
-min{(| + 1)2, 4(p - 2),2(p + n)}. Since A™"VAoi < 2, A02/A01 < 2 and 
A20/A01 < 2 are easily checked in the relevant ranges p G ]2, 6], p G ]6, n + 4] 
and p > n + A, Corollary 11.51 yields an asymptotic expansion weighted by 
the denominator r~'^(p(r, 0)/p(r, ^/))'^+2/("+p) where 7 = A/Aqi is given by 

(dS]), and 5 = ^(f - 1 - v/(| + l)2 + A) agrees with ([TT]). The sum in this 
asymptotic expansion consists of only one term, and it corresponds to the 
eigenvalue Aqi- Thus 

pjr, y) ^ ^ u,,{{l + 2pr)-Py) ( p,{r,y) \ '^^"^''^ 



PB{r,y) ^ B{l + 2pr) \p,{r,Q) ) ^ ' 

as r — 00, where uu{x) / [B + jicp) belongs to the one-dimensional eigenspace 
of Aqi, and the error is measured in an L°° norm with the desired weight. 
Similarly Corollary 11.51 asserts 

+ 0(l/r^). 



Pb{t -ro,y) _ , Mii((l + 2pr) ^y) / pg(r,y) \ 



in the same norm, since the eigenspace is one-dimensional. If tq can be chosen 
to make c(ro) = 1, then subtracting these identities and using the fact ([HD that 
||pB(r — To, y)/ps(r, 1/) — l||Ltx>(j^n) — > will conclude the proof of Corollary II. 71 

On the other hand, from definition ([7]) we compute 

Ps(r-ro,0) ( \ ^ 2p[T - tq)\'^^ ^ , r i ^ (m \ 

= 1 + n/3ro/r + o(l/r) , 



Ps(r,0) V l + 2pr 

from which we deduce that c{tq) does not generally vanish but depends linearly 
on tq. Thus a suitable choice of tq yields c{tq) = 1 to complete the proof. □ 

Proof of Corollary \1.1 (A In the range m > m„+4 and n > 2, the next largest 
spectral value after A20 is An (and not Aq°'^*). Taking A = An and observing 
that the spectral ratio An/Aoi = 2 — {p + A — n)/2p is strictly less than 2 in 
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the range p > n + 4 oi interest, as in the preceding proof we can choose Tq to 
obtain 



2 



P^{t, y) (1 + 2pr)(3p+"-4)/2p \p^{^r, 0) 

as r — 7- oo, where the error is measured in an norm with the desired weight 
according to Theorem II. 2 [ and U2o{y) is a homogeneous harmonic polynomial 
of degree 2 as in Corollary 111.31 

A number of authors starting with Titov and Ustinov |12] and Tartar (circa 
1986, unpublished but cited in [45j) and including two of us [21], have indepen- 
dently observed that the family of functions UB('E~^^'^y) det E~^/^ parametrized 
by positive definite symmetric matrices S > form a invariant manifold of 
dimension n{n + l)/2 under the porous medium and fast diffusion dynamics 
(dJ). Each solution p{T,y) = UB{Jl~^/'^{T)y) detS~^/^(r) satisfies 

p{r - fo, y) . , M2o((l + 2pTy^y) f pjj^^ ^ , ^ f-y^ 



Psir, y) (1 + 2pr)(3p+"-4)/2p ^p^(r, 0) 

as above. Since S(r) is proportional to the moment of inertia tensor y ® 
yp{T, y)dy, Corollary 111.31 shows the traceless part Eq of E(r) is indepen- 
dent of r and can be selected to make -1120 = ^20- The evolution equa- 
tion (dc7/(ir)P+" = cb det S(r) for E(r) = Eq + a(r)/ > is from Corol- 
lary 4 of [2IJ. Subtracting the two equations above and using the fact that 
||p(r — To, ?/)/Ps(t, t/) — 1||l°=(R") — )■ as in ([8]) yields the desired limit f|T8|) . 
after translating the solution p in time by tq — tq. □ 



12. Appendix: Pedestrian derivation of all Schauder Estimates 

Although the announced results have now been established, for the sake of 
convenience, we include this appendix giving a self-contained derivation of all 
Schauder estimates for linear equations that we have relied on. 



Self-contained proof of Lemma l6J[ By superposition, the estimates 
sembled from four cases, in each of which exactly one of the quantities /, b, c, vo 
is non-zero. 

We repeat and modify the arguments from [34j, IV, appropriately to match our 
situation, denoting by T the heat kernel: T(t,x) := (47rt)~"/^ exp[— |a;p/4t]. 

We estimate the spatial Holder quotients for the contribution from /: 
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In this estimate, B{2r,x) denotes the ball centered at x with radius 2r = 
2\x — x'\. We have 

v{t,x)-v{t,x')= f f d^jr{t-T,x-y){P{T,y)-r{T,x))dydT 

Jo JR" 

- /* / ajr(t -T,x'- y) (r (r, y) - r{r, x')) dy dr 

Jo JR" 

= / / dlr(t-r,x-y){r{r,y)-r{r,x))dydr 

Jo JB{2r,x) 

- f I dlV{t -T,x'- y) (r (r, y) - P{t, x')) dy dr 

Jo J B{2r,x) 

+ [ [ (ajr(i -r,x-y)- dlr{t -r,x'- y)) (r, y) - ^ (r, x)) dy dr 

Jo J B{2r,xY 

- f {P{t, x) - P{t, x')) f ajr(t -T,x'- y) dy dr . 

Jo J B{2r,xY 

We note: If |x — y| < 2r, then \x' — y\< 3r. In the third term, if |a; — |/| > 2r, 

then \ax + (1 — cf)x' — y\> \\x — y\ for any a e [0, 1]. Therefore 

(114) 



\v{t,x) -v{t,x')\ < 




dfT{t -r,x-y)\\x- y^ [P{r)l;a dy dr 



J B{2r,x) 



+ 
+ 



J B{2r,x) 



\d%T(t -T,x'-y)\ \x' - si" [J"(t)Uo. dy dr 



^0 JB(2r,xY 



r,x-y)- dfT{t -r,x'-y)\\x- yl" [r{r)]x;a dy dr 



'B{2r,xY 

+ I \x-x'r[P{r)l.,, 



dlV{t-r,x' -y)dy 



dr 



< c 
+ C 
+ C 




J B{2r,x) 



'B{2r,xY 

{t - r)-t-^ exp [-C^] \x-yr M].;a dy dr 




'0 J B(3r,x') 



{t - r)-t-i exp[-C(^] \x' - yr [P {r)Ua dy dr 

- x'\{t - r)-t-i exp[-C^] \x - yr [r{r)Uadydr 
+ C /V-^r[rWka / \Vr{t-r,x'-y)\dS{y)dr. 

Jo JdB{2r.x) 




J B{2r,xY 



'dB{2r,x) 

Now if we estimate [f{r)]x-a < ||/||c«(R5,) and in the first three integrals eval- 
uate the time integration first, using that, for /c > 1, 
(115) 

/ (t-r)-^exp[-^] dr = {A^f-'' / s-'^expf-l] ds < C{A^f-\ 

Jo Jo 
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we conclude that each term comes to the same estimate, namely we get 

\v{t,x)-v{t,x')\ < C||/||c:<.(Rj) Ix-x'r • 

In comparison, if we estimate [f{T)]x;a < ''"~"''^II/IIc"(R")' obtain the fol- 
lowing time integral instead, and we split it in the middle, obtaining 

't/2/ 

pt/2 

(116) -a6[t/2,t]^ ' Jo 

+ (V2)-"/'(A')'-*^ / s-'' exp[-l/s] ds 
Jo 

with the same estimate for each summand. This results in 

\vit,x)-vit,x')\ < C||/||^„(i,„) \x-x'r . 

We now estimate the time Holder quotients for the contribution from /: 
We assume t' > t and let t' — t —: d. Then 

v{t',x)-v{t,x)= f [ df^T{t'-T,x-y){riT,y)-r{T,x))dydr 
Jo Jh" 

^ f f d^^T{t-T,x-y){P{T,y)-P{T,x))dydT 
Jo JR" 

^ f [ dlr{t' -T,x-y) (r (r, y) - ^ (r, x)) dy dr 

J{t-d)+ Jb." 

- f [ df^r{t -T,x-y) (r (r, y) - ^ (r, x)) dy dr 

J{t-d)+ JR" 
r{t-d)+ r 

+ / / [dint' -T,x-y)- dint - r, X - y)] {P{t, y) - ^ir, x)) dy dr . 
Jo -/R" 



Hence 
(117) 

\v{t',x) -v{t,x)\ < 
ft' 



C f f {f- T)-t-^ exp [-C^] \x-yr dy [/(t)],;„ dr 

J{t-d)+ Jw 

+ C t [ {t- r)-t-i exp [-C^] \x-yr dy [/(r)],;„ dr 

J{t~d)+ 7r" 
f{t-d)+ p 

+ C (t'-t) / (r-T)-i-2exp[-C^]|a;-yrdy[/(T)Udr, 

Jo JR" 



where t* G [t,t']; and therefore, in the last domain of integration, we have 
— r) < t^: — T < t' — T. Evaluating the space integrals directly first, we 
obtain 

\v{t' , x) — v(t, x)\ < 

ft' 

(TTo\ C {t' — T)2 -"^ [/(r)]^,;^ dr + same with t instead of t' 

+ C r \t'-t){t'-r)^-'[f{r)l,^dT. 
Jo 



Each term is dominated by sup, ||/(r) ||c;"(R") < C{t' -t)''/^\\f\\c^^n-), 

as desired. 

If we use the weighted norms in psjl instead, we obtain 



\v{t',x) -v{t,x)\ < 

rt' 

C 



J{t-d)+ 

;il9) +C f (t-r)t-V--/2||/||^.(^„)rfr 



{t-d)+ 

+ C{t'-t) '^"(t-r)t-V--/2||/||^.(^„)dr 



< C{t' - t)"/2t-«/2 



To justify the last estimate for the first two integrals, we distinguish two cases: 
If > it, our estimate = ((i/t)"/^ is weaker than a constant, whereas 

the integrals can be extended to a lower limit r = and readily estimated by 
a constant. If d < |t, the term 

— (|^) °^^) 3.nd the other factor can 
be integrated and found to be bounded by 0{{t' — t)"/^). 

For the third integral, in the case that d > we argue that (t — r)t^2 ^ 
— (^;f'_^^f-2^ ^Yie coefficient of ||/||* from the third term is dominated 

by C{t' - t)^-'^{t - c/)+ ^ This quantity is again bounded by a constant, 
whereas our claimed estimate is weaker than a constant. On the other hand, if 
d < |t, we split the integral in the middle, estimate the bounded factor under 
each integral by its maximum (at r = |t for each integral) and integrate the 
remammg factor. Now d fj'^ ...dT< Cd/t < C{d/tYl'^ and d ...dT< 
C{d/t)°'/2, hence the desired estimate. 

We now estimate the supremum norm for the contribution from /: 
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This is the easy estimate 

Ht,x)\< f ! \dlT{t~T,x-y)\\r{T,y)-r{T,x)\dydT 

Jo JR" 

(120) <C f [ {t-r)-'^-'exp[-C\^]\x-yrdy[f{T)l,^dT 

Jo 

< C [\t - r)t-i[/(r)].;, dr < Cr/'\\f\\c.^n-) ■ 



Since t is bounded, we have estimated \v{t, x) \ by the Holder norm of /. In the 
same way, we can also estimate t°'^'^\v{t, x)\ < C'^"^^||/|Ic'q(rj) — ^II/IIc"(R5.)- 

We now estimate the spatial Holder quotient contributed from b: 

They can be estimated as in ( I114p . but without splitting the space integral, as 

\v{t,x)-v{t,x')\ ^ I' I |5^r(t-r,y)|[6V)Wrf2/rf^ 

\X X \ Jq Jjfn 

^^^^^ <C [\t-r)-^b{T)l,^dT 

Jo 

< min{Ct^||6||c<^(Rn), Ct^-f } . 

We now estimate the time Holder quotients contributed from b: 

For the time Holder quotients, we can argue as in (11181) . only with an extra 
power of (t' — r)^/^ under the integrals. We get 

|f(t',x) -v{t,x)\ < C((t' - t)t+i + (t' _ t)(t' _ t)f-3) sup[6(r)]^.„ 

T 

<C(t'-t)tT^||6||c..(R.). 

If we use the weighted norms || • ||* instead, we get by modification of (I119p : 

\v(t', x) — v{t, x)\ < 

C^II&llc^(Rf) r it' - rr/'-^'r'-/' dr 
J{t-d)+ 

^'"'^ +C\\b\\h.,n-J' {t-rr/^-y^r-/^dr 

J{t-d)+ 

+ C(t'-t)||6|rc..(H.) r '^\t-r)^^-lT-'^dT. 
Jo 

The same splitting argument as for (11191) proves that each term is dominated 
by C{t' - t)"/^+^/^t-/2||6|r^.(^„) < CT^^t' - tr'H-l^b\\},.^^.^y 

We now estimate the supremum norm contributed from b: 
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As in fll20p . we conclude \v{t,x)\ < CU^^^ \\b\\c^(Rj^) and \v{t,x)\ < CU ||&||^«(r5;)- 
An even simpler version of the same estimate yields |f()f:,x)| < Ct5 ||6||^oo(pt^) 
for later use in proving (jH]). 

We now estimate the space and time Holder quotients contributed from c: 

In contrast to b and /, we cannot follow the paradigms of ( 1114^ and ( ]117p 
here, because F carries no space derivative here. This impacts the results of 
['y]t;a/2 since we cannot benefit from using spatial Holder quotients of c. We 
get 



\v{t,x) - v{t,x')\ ^ 



/ / r(t - r,?/) [c(r)]^.„rf?/dr < t||c||c-{R^) . 
Jo Jh" 



Similarly, W^[f (t)]^;^ < ^l|c||^<»(R,n). 

For the time Holder quotients, we have (with t' > t) 



R" 



c(t' -T,x-y)-c{t-T,x-y)\ 



+ / / r(r,y) —-—75 dy dr 



jh' 



t' - 



) • 



With |v(t,a;)| < ]j;a/2, the c-estimate in f H5]) follows immediately. The 
example c = 1, v{t,x) = t shows that the estimate is optimal. 

Estimating the time Holder quotients in terms of the weighted norm, we get 
analogously 

\v{t',x) -V{t,x)\ 1 /■*' II s|| , 
1^/ _ t\»/2 ^ _ ^|a/2 -r)|Uoo(R„)rfr 







<C{t — "||c||cQ(Rn) + Ct^ °^^l|c||c«(R5,) • 



Hence i"/2 Mij£l_^(i£)i < qj'^ °''^||c||^„j-j^„-); again, as a corollary, we get 
\v{t',x) \ < CT\\c\\*^a(^^ny Thls glvcs the c estimate in (HSj) . 

The estimates for the contributions from t;o have been given in Sec. [6l Equa- 
tions (117D,(llH]),(linD aheady. 

Likewise, the estimates for the 

contributions of ||&||l°° and ||c||ioo to the alternate bound (jH]) have been given 
in Sec. |6] already. 

This concludes the proof of Lemma 16.11 □ 
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